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ABSTRACT 

We present the results of the Gemini Deep Planet Survey, a near-infrared adaptive optics search for giant 
planets and brown dwarfs around nearby young stars. The observations were obtained with the Altair adaptive 
optics system at the Gemini North telescope and angular differential imaging was used to suppress the speckle 
noise of the central star. Detection limits for the 85 stars observed are presented, along with a list of all faint 
point sources detected around them. Typically, the observations are sensitive to angular separations beyond 
0.5" with 5(T contrast sensitivities in magnitude difference at 1.6 //m of 9.5 at 0.5", 12.9 at 1", 15.0 at 2", 
and 16.5 at 5". For the typical target of the survey, a 100 Myr old KO star located 22 pc from the Sun, the 
observations are sensitive enough to detect planets more massive than 2 Mjup with a projected separation in the 
range 40-200 AU. Depending on the age, spectral type, and distance of the target stars, the detection limit can 
be as low as ^^1 Mjup. Second epoch observations of 48 stars with candidates (out of 54) have confirmed that all 
candidates are unrelated background stars. A detailed statistical analysis of the survey results, yielding upper 
limits on the fractions of stars with giant planet or low mass brown dwarf companions, is presented. Assuming 
a planet mass distribution dn/dm cx m"'-^ and a semi-major axis distribution d«/dfl cx a~\ the 95% credible 
upper limits on the fraction of stars with at least one planet of mass 0.5-13 Mjup are 0.28 for the range 10- 
25 AU, 0.13 for 25-50 AU, and 0.093 for 50-250 AU; this result is weakly dependent on the semi-major axis 
distribution power-law index. The 95% credible interval for the fraction of stars with at least one brown dwarf 
companion having a semi-major axis in the range 25-250 AU is 0.019;!;j]'Q[5, irrespective of any assumption on 
the mass and semi-major axis distributions. The observations made as part of this survey have resolved the 
stars HD 14802, HD 166181, and HD 213845 into binaries for the first time. 

Subject headings: Planetary systems — stains: imaging — binaries: close — stains: low-mass, brown dwaifs 



1. INTRODUCTION 

More than 200 exoplanets have been discovered over the 
last decade through precise measurements of variations of the 
radial velocity (RV) of their primary star. Besides establish- 
ing that at least 6-7% of FGK stars have at least one giant 
planet with a semi-major axis smaller than ^5 AU (Marcy 
et al. 2005), the profusion of data following from the RV dis- 
coveries has propelled the field of giant planet formation and 
evolution into an unprecedented state of activity. For a re- 
view of the main characteristics of the RV exoplanets, the 
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reader is referred to Udry et al. (2007); Butler et al. (2006); 
Marcy et al. (2005). Besides the RV technique, the photo- 
metric transit method has lead successfully to the discovery 
of new exoplanets on small orbits (e.g. Konacki et al. 2003; 
Alonso et al. 2004; Cameron et al. 2007) and has provided the 
first measurements of the radius and mean density of giant ex- 
oplanets (e.g. Charbonneau et al. 2000). Very recently, a few 
exoplanets have been detected by gravitational microlensing 
(Bond et al. 2004; Udalski et al. 2005; BeauUeu et al. 2006; 
Gould et al. 2006); these planets have separations of ~2- 
5 AU. Notwithstanding their great success in finding planets 
on small orbits, these techniques cannot be used to search for 
and characterize planets on orbits larger than ~10 AU. As a 
result, the population of exoplanets on large orbits is currently 
unconstrained. 

The two main models of giant planet formation are core 
accretion (Pollack et al. 1996) and gravitational instability 
(Boss 1997, 2001). In the core accretion model, solid particles 
within a proto-planetary disk collide and grow into solid cores 
which, if they become massive enough before the gas disk dis- 
sipates, trigger runaway gas accretion and become giant plan- 
ets. Models predict that the timescale for formation of a planet 
like Jupiter through this process is about 5 Myr (Pollack et al. 
1996), or about 1 Myr if migration of the core through the disk 
is allowed as the planet forms (Alibert et al. 2005). These 
timescales are comparable to or below the estimated proto- 
planetary dust disk lifetime (~6 Myr, Haisch et al. 2001) and 
gas disk lifetime (<10 Myr, Jayawardhana et al. 2006). For- 
mation through core accretion is strongly dependent on the 
surface density of solid material (hence [Fe/H]) in the proto- 
planetary disk , precluding formation of Jupiter mass plan- 



2 



Lafreniere et al. 



ets at distances greater than 15-20 AU (e.g. Pollack et al. 
1996; Ida & Lin 2004), where the low density of planetesi- 
mals would lead to prohibitively long formation timescales. 
Neptune mass planets can be formed out to slightly larger dis- 
tances and can further migrate outward owing to interaction 
with the disk. 

In the gravitational instability model, small instabilities in a 
proto-planetary disk grow rapidly into regions of higher den- 
sity that subsequently evolve into spiral arms owing to Ke- 
plerian rotation. Further interactions between these spiral 
arms lead to the formation of hot spots which then collapse 
to form giant planets. The range of orbital separation over 
which this mechanism may operate efficiently is not yet clear. 
Some studies indicate that it may lead to planet formation only 
at separations exceeding ~100 AU (Whitworth & Stamatel- 
los 2006; Matzner & Levin 2005), where the radiative cool- 
ing timescale is sufficiently short compared to the dynamical 
timescale, while others have been able to produce planets only 
at separations below 20-30 AU (Boss 2000, 2003, 2006). 

A few other models are capable of fornaing giant planets on 
large orbits directly. One such mechanism is shock-induced 
formation following collision between disks (Shen & Wad- 
sley 2006). In this model, the violent collision of two proto- 
planetary disks triggers instabilities that lead to the collapse of 
planetary or brown dwarf (BD) mass clumps. Results of nu- 
merical simulations indicate that planets and BDs may form 
at separations of several tens of AU or more through this pro- 
cess (Shen & Wadsley 2006). The competitive accretion and 
ejection mechanism that was proposed initially to explain the 
formation of BDs (Reipurth & Clarke 2001) could also form 
planetary mass companions on large orbits, as suggested by 
the results of recent simulations by Bate & Bonnell (2005). 

Even in a scenario in which all giant planets form on small 
orbits, through either core accretion or gravitational collapse, 
a significant fraction of planets could be found on stable or- 
bits of tens of AU because of outward orbital migration. In- 
deed, numerical simulations have shown that gravitational in- 
teractions between planets in a multi-planet system may send 
one of the planets, usually the least massive one, out to an 
eccentric orbit of semi-major axis of tens to hundreds of AU 
(Chatterjee et al. 2007; Veras & Armitage 2004; Rasio & Ford 
1996; Weidenschilling & Marzari 1996). This process could 
be involved frequently in the shaping of the orbital parame- 
ters of planetary systems as we have learned from RV surveys 
that multi-planet systems are common, representing ~14% of 
known planetary systems (Marcy et al. 2005). Outward mi- 
gration of massive planets can be induced also by interactions 
between the planet and the gaseous disk; the simulations of 
Veras & Armitage (2004) reveal that this process is capable of 
carrying Jupiter mass planets out to several tens of AU. Sim- 
ilarly, angular momentum exchange between two planets (or 
more), achieved through viscous interactions with the disk, 
could drive the outer planet to a separation of hundreds of AU 
(Martin et al. 2007). Outward planet migration can result fur- 
ther from interaction of the planet with the soUd particles in 
the disk after the gas has dissipated (e.g. Levison et al. 2007); 
there is in fact strong evidence that this mechanism has played 
an important role in the Solar system (Fernandez & Ip 1984; 
Malhotra 1995; Hahn & Malhotra 2005). Based on numerical 
simulations, it is Ukely that all giant planets of the Solar sys- 
tem formed interior to -^15 AU and migrated outward (except 
Jupiter) to their current location (Tsiganis et al. 2005). 

From an observational point of view, there is some evidence 
that planets on large orbits may exist. Many observations of 



dusty disks around young stars, made either in emitted light 
(e.g. Vega, e Eri, Fomalhaut; Holland et al. 1998; Greaves 
et al. 1998) or in scattered fight (e.g. HD 141569, HR 4796, 
Fomalhaut; Augereau et al. 1999; Weinberger et al. 1999; 
Schneider et al. 1999; Kalas et al. 2005), have unveiled asym- 
metric or ring-like dust distributions. These peculiar mor- 
phologies could arise from gravitational dust confinement im- 
posed by one or more (unseen) giant planets on orbits of tens 
to hundreds of AU. In fact, detailed numerical simulations of 
the effect of giant planets on the dynamical evolution of dusty 
disks have been able to reproduce the observed morphologies 
with remarkable agreement (Ozernoy et al. 2000; Wilner et al. 
2002; Deller & Maddison 2005). Typically, Jupiter mass plan- 
ets on orbits of ~60 AU are needed to reproduce the observa- 
tions, although in some cases less massive planets (similar to 
Neptune) may be able to reproduce the observed features. 

In the last few years, there have been a few discoveries 
of planetary mass or low-mass BD companions located be- 
yond several tens of AU, in projection, from their primary: an 
~8Mjup companion 40 AU from the BD 2M 1207-3932(Mo- 
hanty et al. 2007; Chauvin et al. 2005a, 2004), a ~25 Mjup 
companion 100 AU from the T Tauri star GQ Lup (Marois 
et al. 2007a; Seifahrt et al. 2007; Neuhauser et al. 2005), a 
~12 Mjup companion 210 AU from the young star CHXR 73 
(Luhman et al. 2006), a ~25 Mjup companion 240 AU from 
the young BD 2M 1101-7732 (Luhman 2004), a ~12 Mjup 
companion 260 AU from the young star AB Pic (Mohanty 
et al. 2007; Chauvin et al. 2005b), a 7-19 Mjup companion 
240-300 AU from the young BD Oph 1622-2405 (Luhman 
et al. 2007a; Close et al. 2007; Jayawardhana et al. 2006), an 
~11 Mjup companion 330 AU from the T Tauri star DH Tau 
(Luhman et al. 2006; Itoh et al. 2005), and a --21 Mjup com- 
panion 790 AU from the star HN Peg (Luhman et al. 2007b). 
These discoveries might indicate that more similar compan- 
ions, and less massive ones, do exist and remain to be found. 

Perhaps even more compelling is the fact that the number 
of exoplanets found by RV surveys increases as a function of 
semi-major axis for the range 0.1-3 AU (Butler et al. 2006); 
these surveys are incomplete at larger separations. Conserva- 
tive extrapolation suggests that there may be at least as many 
planets beyond 3 AU as there are within (Butler et al. 2006). 
In fact, long-term trends in RV data have been detected for 
about 5% of the stars surveyed (Marcy et al. 2005), suggest- 
ing the presence of planets between 5 AU and 20 AU around 
them. 

Given all of the considerations above, it is clear that a de- 
termination of the frequency of giant planets as a function of 
orbital separation out to hundreds of AU is necessary to elu- 
cidate the relative importance of the various modes of planet 
formation and migration. Direct imaging is currently the only 
viable technique to probe for planets on large separations and 
achieve this goal. However, detecting giant planets directly 
through imaging is very difficult due to the angular proxim- 
ity of the star and the very large luminosity ratios involved. 
Currently, the main technical difficulty when trying to image 
giant planets directly does not come from diffraction of light 
by the telescope aperture, from fight scattering due to resid- 
ual atmospheric wavefront errors after adaptive optics (AO) 
correction, nor from photon noise of the stellar point spread 
function (PSF), but rather from light scattering by optical im- 
perfections of the telescope and camera that produce bright 
quasi-static speckles in the PSF of the central star. These 
speckles are usually much brighter than the planets sought 
after. More in depth discussions of this problem, as well as 
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Fig. 1 . — Distribution of distance, spectral type, and age of tiie target stars. For tiie age distribution, each star was distributed over all the age bins according to 
the fraction of their estimated age interval faffing inside each bin. 



possible venues to circumvent it using current instrumenta- 
tion, can be found in Lafreniere et al. (2007); Hinkley et al. 
(2007); Marois et al. (2006, 2005); Masciadri et al. (2005); 
Biller et al. (2004); Schneider & Silverstone (2003); Marois 
et al. (2003); Sparks & Ford (2002); Marois et al. (2000); 
Racine et al. (1999). As AO systems continue to improve 
and eventually achieve Strehl ratios above ^90%, the light 
diffracted by the telescope aperture will become more impor- 
tant compared to scattered light and the use of a coronagraph 
will be mandatory. But even then, after removal of diffracted 
Ught by the coronagraph, high-contrast imaging appUcations 
will likely be limited by residual quasi-static speckles (e.g. 
Macintosh et al. 2006). 

Many direct imaging searches for planetary or brown dwarf 
companions to stars have been done during the last five years, 
see for example Biller et al. (2007); Chauvin et al. (2006); 
Metchev (2006); Lowrance et al. (2005); Masciadri et al. 
(2005); McCarthy & Zuckerman (2004); and Luhman & 
Jayawardhana (2002) for searches carried out in J, H, or K, 
or Kasper et al. (2007) and Heinze et al. (2006) for searches 
made in L' or M' . Depending on the observing strategy em- 
ployed, the properties of the target stars, and the characteris- 
tics of the instrument used, each of these surveys was sensitive 
to a different regime of companion masses and separations. 
Typically, these surveys have reached detection contrasts of 
10-13 mag for angular separations beyond l"-2", sufficient 
to detect planets more massive than ^5 Mjup for targets aged 
~100 Myr. Unfortunately, rigorous statistical analyses allow- 
ing derivation of clear constraints on the population of planets 
in the regimes of mass and separation to which these surveys 
were sensitive are orily begirming to be reported in the Utera- 
ture^ ; an assessment of the current status of knowledge is thus 
rather difficult to make. Nonetheless, it is fair to say that the 
population of planets less massive than ^5 Mjup, having orbits 
with a semi-major axis of tens to hundreds of AU, is poorly 
constrained. 

In this paper we report the results of the Gemini Deep 
Planet Survey (GDPS), a direct imaging survey of 85 nearby 
young stars aimed at constraining the population of Jupiter 
mass planets with orbits of semi-major axis in the range 10- 
300 AU. The selection of the GDPS target sample is explained 
in §2, and the observations and data reduction are detailed in 
§3. The detection Umits achieved for each target are then pre- 
sented in §4 along with all candidate companions detected. A 

' In addition to the present work, analyses by Nielsen et al. (2007) and 
Kasper et al. (2007) have become available during the review process of this 
manuscript. 



Statistical analysis of the results allowing determination of the 

maximum fraction of stars that could bear planetary compan- 
ions is presented in §5. Concluding remarks follow in §6. 

2. TARGET SAMPLE 

In light of the luminosity ratio and angular separation prob- 
lem highlighted above, the list of target stars was assembled 
mainly on the basis of young age and proximity to the Sun, 
the latter yielding a larger angular separation for a given phys- 
ical distance between the star and an eventual planet. Equiv- 
alently, a given detection threshold is achieved at a smaller 
physical separation for a star closer to the Sun, and planets on 
smaller orbits can be detected. Additionally, for angular sep- 
arations where planet detection is limited by sky background 
noise or read noise, lower mass planets can be detected around 
a star closer to the Sun as their apparent brightness would be 
larger Giant planets are intrinsically more luminous at young 
ages and fade with time (e.g. Marley et al. 2007; Baraffe et al. 
2003; Burrows et al. 1997); therefore, for a given detection 
threshold, observations of younger stars are sensitive to plan- 
ets having a lower mass. The proximity and age criteria used 
in building the target list thus maximize the range of mass and 
separation over which the survey is sensitive. 

The target stars were selected from three sources: (1) Ta- 
bles 3 and 4 of Wichmann et al. (2003), which list nearby 
stars with an estimated age below or comparable to that of 
the Pleiades (^100 Myr), based on measurements of lithium 
abundance, space velocity, and X-ray activity; (2) Tables 2 
and 5 of Zuckerman & Song (2004b), which list members of 
the /3 Pictoris (~12 Myr) and AB Doradus (~50 Myr) mov- 
ing groups respectively; and (3) Tables 2 and 5 of Montes 
et al. (2001b), which list late-type single stars that are possible 
members of the Local Association (Pleiades moving group, 
20-150 Myr) and IC 2391 supercluster (35-55 Myr) respec- 
tively, based on space velocity measurements. The stars Usted 
in Montes et al. (2001b) were initially selected based on vari- 
ous criteria indicative of youth, such as kinematic properties, 
rotation rate, chromospheric activity, lithium abundance, or 
X-ray emission, but for many of these stars the space velocity 
is the only indication of youth as other measurements are ei- 
ther unavailable or inconclusive; the young age of such stars 
is therefore uncertain. This uncertainty will be taken into ac- 
count in our statistical analysis (§5). A few stars known to 
have a circumsteUar disk were added to these lists. 
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GDPS TARGET SAMPLE 







Names 




(X 


S 


Spectral 


H 


Dist." 


fi(X COS 5^ 




[Fe/H] 


[Fe/H] 


Age 


Age 


Notes 


HD 


GJ 


HIP 


Other 


(J2000) 


(J2000) 


Type 


(mag) 


(pc) 


(mas/yr) 


(mas/yr) 




Ref. 


(Myr) 


Ref. 




166 


5 


544 


- 


00i"06"36?78 


+29°01'17"4 


KO V 


4.63 


13.7 


379.9 


-178.3 


0.18 


V05 


150-300 


G98,G00,F04,L06 


Her-Lyr 


691 


- 


919 


V344 And 


00i"ll"22?44 


+30°26'58"5 


KO V 


6.26 


34.1 


210.7 


35.3 


0.32 


V05 


50-280 


W03.W04.C05 


- 


1405 


- 


- 


PWAnd 


00'"18""20?90 


+30°57'22"0 


K2 V 


6.51 


30.6'' 


143.7'' 


-171.5'' 


-0.78 


N04 


50-50 


Z()4b 


AB Dor 


5996 


- 


4907 


- 


0l''02"'57?22 


+69°13'37"4 


G5 V 


5.98 


25.8 


223.9 


-148.4 


-0.28 


K02 


100-650 


M01a,G03 


?LA 


9540 


59A 


7235 


- 


Ol'"33"15?81 


-24°10'40"7 


KO V 


5.27 


19.5 


271.9 


-159.5 


-0.02 


V05 


100-1350 


M01a,W04 


?LA 


10008 


- 


7576 


- 


0l'"37"35?47 


-06°45'37"5 


G5 V 


5.90 


23.6 


171.0 


-97.7 


- 


- 


150-300 


F04,L06 


Her-Lyr 


- 


82 


9291 


V596 Cas 


0l'"59"23?51 


+58°31'16"1 


dM4e 


7.22 


12.0 


321.8 


-195.3 


- 


- 


50-50 


M06 


a Per 


14802 


97 


11072 


kap For 


02h22"32?55 


-23°48'58"8 


G2 V 


3.71 


21.9 


197.3 


-4.4 


-0.04 


T05b 


5000-6700 


L99,W04,B99 


m 


16765 


- 


12530 


84 Get 


02'"41"'14?00 


-00°41'44"4 


F7 V 


4.64 


21.6 


218.9 


-125.3 


-0.27 


N04 


100^00 


M01a,H98,F95 


LA,m 


17190 


112 


12926 


- 


02''46"'15?21 


+25°38'59"6 


Kl V 


6.00 


25.7 


238.7 


-149.9 


-0.11 


V05 


50-3500 


M01a,W04 


?1C 2391 


17382 


113 


13081 


- 


02'>48"'09?14 


+27°04'07;'l 


Kl V 


5.69 


22.4 


264.2 


-127.8 


0.12 


T05b 


50-600 


M01a,W04 


?1G 2391 


17925 


117 


13402 


EPEri 


02''52'"32?13 


-12°46'li;'0 


K2V 


4.23 


10.4 


398.1 


-189.6 


0.18 


V05 


40-128 


W03,M01a,M01b,L99 


LA 


18803 


120.2 


14150 


51 Ari 


03''02'"26?03 


+26°36'33;'3 


G8V 


5.02 


21.2 


234.2 


-168.2 


0.11 


V05 


800-3600 


W04,T05a 


- 


19994 


128 


14954 


94 Get 


03''12'"46?44 


-01°11'46;'0 


F8 V 


3.77 


22.4 


193.4 


-69.2 


0.19 


V05 


800-3500 


W04,T05a 


m 


20367 


- 


15323 


- 


03''17'"40?05 


+31°07'37;'4 


GOV 


5.12 


27.1 


-103.1 


-56.6 


0.17 


S04a 


50-150 


W03 


- 


- 


- 


- 


2E759 


03''20"'49?50 


-19°16'10f'0 


K7V 


7.66 


27.0 


90.8 


^3.8 


- 


- 


50-150 


M01a,L06 


LA 


22049 


144 


16537 


eps Eri 


03''32"'55?84 


-09°27'29;'7 


K2V 


1.88 


3.2 


-976.4 


18.0 


-0.03 


V05 


530-930 


SOO 


- 


- 


- 


17695 


- 


03''47'"23?35 


-01°58'19;'9 


M3e 


7.17 


16.3 


186.7 


-271.8 


- 


- 


80-120 


L06 


LA (B4) 


25457 


159 


18859 


- 


04h02'"36?74 


-00°16'08"1 


F6 V 


4.34 


19.2 


151.2 


-252.0 


0.02 


T05b 


80-120 


L06 


LA (B4) 


283750 


171.2 


21482 


V833 Tau 


04h36™48?24 


+27°07'55"9 


K2 V 


5.40 


17.9 


232.4 


-147.1 


- 


- 


50-150 


W()3 


- 


30652 


178 


22449 


1 Ori 


04'"49""50?41 


+06°57'40"6 


F6 V 


1.76 


8.0 


463.4 


11.6 


0.03 


V05 


50-500 


M01a,H99 


IC 2391 


- 


182 


23200 


V1005 Ori 


04'>59"34?83 


+01°47'00"7 


Ml V 


6.45 


26.7 


37.1 


-93.9 


- 


- 


10-50 


F98,B99 


- 


- 


234A 


30920 


V577 Mon 


06''29""23?40 


-02°48'50"3 


M4 


5.75 


4.1 


694.7 


-618.6 


- 


- 


100-3000 


M01a,M03 


?LA,m 


- 


281 


37288 


- 


07'"39""23?04 


+02°11'01"2 


K7 


6.09 


14.9 


-147.9 


-246.6 


- 


- 


150-300 


L06 


Her-Lyr 


- 


285 


37766 


YZCMi 


07'"44"40n7 


+03°33'08"8 


M4.5 V 


6.01 


5.9 


-344.9 


-450.8 


- 


- 


50-50 


Z,M01 


LA 


72905 


311 


42438 


3Uma 


08l>39""ll?70 


+65°01'15"3 


G1.5 V 


4.28 


14.3 


-27.7 


87.9 


-0.09 


T05b 


300-300 


S93,M01b 


UMa 


75332 


- 


43410 


- 


08l"50"32?22 


+33°17'06"2 


F7 Vn 


5.04 


28.7 


-62.2 


-85.0 


0.14 


V05 


50-150 


W03 


- 


77407 


- 


44458 


- 


09l"03"27?08 


+37°50'27"5 


GO 


5.53 


30.1 


-86.3 


-168.8 


0.10 


V05 


10-50 


W03,M01a,M01b 


LA,m 


78141 


- 


- 


- 


09h07"18?08 


+22°52'21"6 


KO V 


5.92 


HA" 


-0.4'' 


-67.6'' 


- 


- 


50-150 


W03 


- 


82558 


355 


46816 


LQHya 


09''32™25?57 


-ll°ll'04f'7 


KO V 


5.60 


18.3 


-248.2 


35.1 


0.33 


V05 


50-100 


W03,M01b 


- 


82443 


354.1 


46843 


DXLeo 


09''32'"43?76 


+26°59'18;'7 


KO V 


5.24 


17.7 


-147.5 


-246.3 


-0.10 


T05b 


50-150 


W03,M01b 


LA 


- 


393 


51317 


- 


10''28'"55?55 


-fO0°50'27;'6 


M2 


5.61 


7.2 


-602.3 


-731.9 


- 


- 


80-120 


L06 


LA(B4) 


90905 


- 


51386 


- 


10''29"'42?23 


+Ol°29'28;'0 


F5 


5.60 


31.6 


-151.4 


-125.3 


0.07 


V05 


50-150 


W03 


- 


91901 


- 


51931 


- 


10''36"'30?79 


-13°50'35f'8 


K2V 


6.64 


31.6 


-164.6 


23.8 


-0.03 


K02 


50-5000 


MOla 


?IC 2391 


92945 


3615 


52462 


- 


10''43"'28?27 


-29°03'5i;'4 


Kl V 


5.77 


21.6 


-215.4 


-48.5 


0.13 


V05 


80-120 


L06,W03,S04b,W04 


LA(B4) 


93528 


- 


52787 


- 


10''47'"31?16 


-22°20'52;'9 


KOV 


6.56 


34.9 


-122.7 


-29.4 


0.11 


K02 


50-150 


W03,S02 


- 


- 


402 


53020 


EE Leo 


10''50'"52?06 


+06°48'29f'3 


M4 


6.71 


5.6 


-804.4 


-809.6 


- 


- 


150-300 


L06 


Her-Lyr 


96064 


- 


54155 


- 


llii04"'41?47 


-04°13'15r'9 


G4 


5.90 


24.6 


-178.0 


-104.1 


-0.01 


T05b 


50-150 


W03 


m 


97334 


417 


54745 


- 


ll''12"'32?35 


+35°48'50;'7 


GOV 


5.02 


21.7 


-248.6 


-151.3 


0.09 


V05 


80-300 


K01,M01a 


?LA 


102195 


- 


57370 


- 


ll''45™42?29 


+02°49'17"3 


KOV 


6.27 


29.0 


-190.3 


-111.4 


0.05 


S06 


100-5000 


MOla 


?LA 


102392 




57494 




ll''47™03?83 


-ll°49'26f'6 


K2 


6.36 


24.6 


-206.3 


-60.7 






100-5000 


MOla 


?LA,m 


105631 


3706 


59280 




12'"09"'37?26 


+40°15'07"4 


KOV 


5.70 


24.3 


-314.0 


-51.3 


0.20 


V05 


1600-1600 


M01b,W04 




107146 




60074 




12''19"'06?50 


+16°32'53f'9 


G2 V 


5.61 


28.5 


-175.7 


-148.3 


-0.03 


V05 


50-100 


W03,Z04a 




108767B 






del Grv B 


12''29"'51?85 


-16°30'55"6 


KOV 


6.37 


26.9 


-210.0 


-139.3 






40-260 


R05,G01,M01a 


LA 


109085 


471.2 


61174 


eta Grv 


12'"32""04?23 


-16°11'45"6 


F2 V 


3.37 


18.2 


-424.4 


-58.4 


-0.05 


N04 


600-1300 


M03,W05 










BD+60 1417 


12'"43"33?28 


+60°00'52"7 


KO 


7.36 


17.7'^ 


-125.2'' 


-66.4'' 






50-150 


W03 




111395 


486.1 


62523 




12''48'"47?05 


+24°50'24f'8 


G5 V 


4.71 


17.2 


-334.6 


-106.1 


0.13 


V05 


600-1200 


H99,W04,T05a 




113449 




63742 




13''03'"49?65 


-05°09'42;'5 


G5V 


5.67 


22.1 


-189.8 


-219.6 


-0.22 


T05b 


80-120 


L06 


LA(B4) 



TAB LEI — Continued 





Names 






s 


Spectral 


H 


Dist." 


cos (^'' 




FFe/Hl 


FFe/Hl 


Aee 


Ase 


Notes 


HI) 


c;.i 


HIP 


Other 


(J2()()()) 


(J2()()()) 


T\ pe 


(mag) 


(pc) 


(niasA r) 


(iiias/yrj 




Ref. 


(Myr) 


Rcf. 




- 


507.1 


65016 


- 


13'"19""40?12 


+33°20'47"5 


M1.5 


6.64 


17.4 


-298.8 


-143.9 


- 


- 


100-5000 


MO la 


?LA 


116956 




65515 


- 


13''25"'45?53 


+56°58'13"8 


G9 V 


5.48 


21.9 


-217.4 


11.2 


0.06 


T05b 


100-500 


M01a,G00 


LA 


118100 


517 


66252 


EQVir 


13''34"'43?21 


-08°20'31"3 


K5 Ve 


6.31 


19.8 


-287.4 


-91.0 


0.00 


C97 


50-50 


L05,M01a 


?LA 


- 


524.1 


67092 


- 


13''45"'05?34 


-04°37'13"2 


K5 


7.33 


25.7 


-159.9 


-95.3 


- 


- 


100-5000 


MOla 


?LA 


124106 


3827 


69357 


- 


14'"11""46M7 


-12°36'42"4 


Kl V 


5.95 


23.1 


-257.5 


-179.5 


-0.10 


V05 


800-1500 


H99,M01a,W04 


- 


125161B 


9474B 




- 


14i,j6m 12^16 


+51°22'34f'7 


Kl 


6.32 


29.8 


-150.0 


89.4 


- 


- 


50-5000 


MOla 


?IC 2391 


129333 


559.1 


71631 


EKDra 


14'"39"00?21 


+64°17'30f'0 


GOV 


6.01 


33.9 


-138.6 


-11.9 


0.16 


V05 


20-120 


M01b,W03,L06 


LA (B4),m 


130004 




72146 


- 


14'>45"'24?18 


+13°50'46f'7 


KO V 


5.67 


19.5 


-232.9 


-225.7 


-0.24 


K02 


100-5000 


MOla 


?LA 


130322 




72339 


- 


14''47'"32?73 


-00°16'53;'3 


KO V 


6.32 


29.8 


-129.6 


-140.8 


0.01 


V05 


770-2400 


W04,S05 


- 


130948 


564 


72567 


- 


14''50"15?81 


+23°54'42;'6 


Gl V 


4.69 


17.9 


144.7 


32.4 


0.05 


V05 


50-150 


W03 


- 


135363 




74045 


- 


15''07"'56?26 


+76°12'02f'7 


G5 


6.33 


29.4 


-131.9 


169.3 


-0.10 


K02 


35-100 


W03,M01a,C05 


IC 2391,m 


139813 




75829 


- 


15''29"'23?59 


+80°27'0i;'0 


G5 


5.56 


21.7 


-217.4 


105.5 


0.14 


V05 


50-150 


W03 


- 


141272 


3917 


77408 


- 


15''48'"09?46 


+01°34'18;'3 


G8V 


5.61 


21.3 


-176.2 


-166.7 


-0.08 


T05b 


150-340 


G98,G00,W04,L06 


LA 


147379B 


617B 


79762 


EWDra 


16''16'"45?31 


+67°15'22;'5 


M3 


6.30 


10.7 


^85.6 


90.7 


0.00 


V04 


100-1000 


MOla, H99 


?LA 


- 


628 


80824 


V2306 Oph 


16''30'"18?06 


-12°39'45y3 


M3.5 


5.37 


4.3 


-93.6 


-1184.9 


-0.25 


COl 


100-5000 


MOla 


LA 


- 




81084 


- 


16''33'"41?61 


-09°33'12;'0 


M0.5 


7.78 


31.9 


-67.4 


-179.7 


- 


- 


80-120 


L06 


LA(B4) 


160934 


4020A 


86346 


- 


17''38"'39?63 


+61°14'16f'l 


K7 


7.00 


24.5 


-31.2 


59.4 


- 


- 


30-50 


Z04c,L06 


AB Dor.m 


162283 


696 


87322 


- 


1 7'" 50™ 34? 03 


-06°03'01"0 


MO 


6.70 


21.9 


-26.1 


-131.4 


- 


- 


100-5000 


MOla, B98 


- 


166181 




88848 


V815 Her 


18''08'"16?03 


+29°41'28f'l 


G6 V 


5.76 


32.3"= 


107.0= 


-31.0= 


-0.70 


N04 


50-150 


W03 


m 


167605 




89005 


LPDra 


18''09'"55?50 


+69°40'49;'8 


K2 V 


6.46 


31.0 


-25.3 


193.9 


0.13 


K02 


50-5000 


MOla 


?IC 2391,m 


187748 




97438 


- 


19h48'"15?45 


+59°25'22;'4 


GO 


5.32 


28.4 


15.8 


116.5 


-0.06 


N04 


50-150 


W03 


- 


- 


791.3 


101262 


- 


20'"31"32?07 


+33 ° 46' 33 "1 


K5 V 


6.64 


26.2 


142.1 


16.6 


- 


- 


50-1000 


M01a,H99 


?IC 2391 


197481 


803 


102409 


AUMic 


20'"45""09?53 


-31°20'27"2 


MO 


4.83 


9.9 


280.4 


-360.1 


- 


- 


8-20 


ZOl 


/3Pic 


201651 




104225 


- 


2l''06"'56?39 


+69°40'28"5 


KO 


6.41 


32.8 


109.0 


66.0 


-0.18 


K02 


50-5900 


M01a,W04 


?IC 2391 


202575 


824 


105038 


- 


2l''16"'32M7 


+09°23'37"8 


K3 V 


5.53 


16.2 


146.5 


-119.1 


0.04 


V05 


100-1000 


M01a,H99 


?LA 


- 


4199 


106231 


LOPeg 


2l'"31"01?71 


+23°20'07"4 


K8 


6.52 


25.1 


134.1 


-144.8 


- 


- 


30-50 


Z04c,L06 


AB Dor 


206860 


836.7 


107350 


HNPeg 


2ll"44™31?33 


+14°46'19"0 


GOV 


4.60 


18.4 


231.1 


-113.5 


-0.02 


V05 


150-300 


L06,G98,G00 


Her-Lyr 


208313 


840 


108156 


- 


2l''54"H5?04 


+32°19'42"9 


KO V 


5.68 


20.3 


210.6 


-233.4 


-0.04 


V05 


100-1000 


M01a,H99 


?LA 


- 






V383 Lac 


22''20'"07?03 


+49° 30' 1 If 8 


Kl V 


6.58 


27.5'' 


93.4'' 


5.0" 


- 


- 


50-150 


W03,M01a,C05 


LA 


213845 


863.2 


111449 


ups Aqr 


22''34n'41?64 


-20°42'29;'6 


F7V 


4.27 


22.7 


221.6 


-146.6 


0.11 


T05b 


150-300 


L06 


Her-Lyr,m 




875.1 


112909 


GTPeg 


22''51'"53?54 


+31°45'15;'2 


M3 


7.13 


14.2 


527.0 


-50.6 






200-300 


L05,M01a 


?IC 2391 




876 


113020 


ILAqr 


22''53'"16?73 


-14°15'49;'3 


M4 


5.35 


4.7 


960.3 


-675.6 






100-5000 


MOla 


?LA 




9809 


114066 




23''06'"04?84 


+63°55'34;'4 


MO 


7.17 


24.9 


171.0 


-58.5 






30-50 


Z04c,L06 


AB Dor 


220140 




115147 


V368 Cep 


23''19'"26?63 


+79°00'12;'7 


Kl V 


5.51 


19.7 


201.3 


71.6 


-0.64 


N04 


50-150 


M01b,W03 


LA,m 


221503 


898A 


116215 




23''32'"49?40 


-16°50'44;'3 


K5 


5.61 


13.9 


343.6 


-217.7 


0.00 


C04 


100-800 


M01a,H99 


?LA 




900 


116384 




23''35'"00?28 


+01°36'19;'5 


K7 


6.28 


19.3 


343.0 


26.8 


-0.10 


COl 


150-250 


Z06 


Ca-Near,m 




907.1 


117410 




23''48'"25?69 


-12°59'14;'8 


K8 


6.49 


27.1 


233.6 


22.1 






150-250 


Z06 


Ca-Near,m 



References. — (Z) B. Zuckcrman, private communicalion.; (B99) Bmado y Navascues el al. 1999; (C05) Carpenter et al. 2005; (COl) Cayrd dc Slrobcl cl al. 2001; (C97) Cajrd dc Slrobcl cl 1997; (C04) Clem cl al. 2004; (F95) Favala cl al. 1995; (F98} Favala cl al. 1998; (,F04} FLiliniiami 201)4, (GOO) Gaidcts cl al. 2000, (Cm) 
Guidos 1998;(G01)Gcrbaldiclal. 2001;(G03) Gray clal. 2003; (H98) Hucnsch et al. 1998; (H99) Hiinsch et al. 1999; (KOI) Kirkpatrick et al. 2001; (K02) Koloncva cl al. 2002; (L99) Lachaume ct al. 1999; (L06) Lopez-Santiago et al. 2006; (L05) Lowranee cl ai. 2005; |M06) Makarov 2006; |M03) Mohanly & Basri 2003; (MOla) Monies 
et al. 2001b; (MOlb) Monies et al. 2001a; (N04) Nordstrom et al. 2004; (R05) Riekeet al. 2005; (SOS) Saffe et al. 2005; (S04a) Santos et al. 2004; (S93) Soderblom & Mayor 1993; (SOO) Song et al. 2000; (S04b) Song el al. 2004; (S06) Sousa ei al. 2006; (T05a) Takeda & Kawanomoio 2005; (T05b) Taylor 2005; (V04) Valdes el al. 2004; 
(V05)Valenti & Fischer 2005; (W03) Wichmann et al. 2003; (W04) Wright et al. 2004; (W05) Wyatt et al. 2005; (ZOl) Zuckcrnian ct al. 2001;(Z04b) Zuckerman & Song 2004b; (Z04a) Zuckerman & Song 2004a; (Z04c) Zuckerman et al. 2004; (Z06) Zuckerman et al. 2006 

Note. — Star is a member of (a Per) a Persei; (AB Dor) AB Doradus; (/3 Pic) (3 Pictoris; (Ca-Near) Carina-Near; (Her-Lyr) Hercules-Lyra; (LA) Local association; (LA (B4)) Local association, subgroup B4; (U Ma) Ursa Major If a question mark precedes the association, the membership is doubtful or based on kinematics only. An 
"m" indicates that the star is a multiple, see § 4.3 for more detail. 

From the 1 lipparcos catalog (Pcrryman & ES A 1997), unless stated otherwise. 
^ From Monli:s i:l a). (2001b). 

From the Tycho catalog (H0g et al. 1997). 
^ From the Tycho-2 catalog (H0g et al. 2000). 
^ Fix)mFekeletal.(2005). 
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From this preliminary compilation, we have retained only 
stars with a distance smaller than 35 pc, and we have excluded 
stars of declination below -32° since observations were to be 
made from the Gemini North observatory. Finally, we have 
further excluded stars indicated to be multiple in Zuckerman 
& Song (2004b). This procedure yielded a list of slightly over 
100 target stars, of which 85 were actually observed. The 
properties of these 85 stars are presented in Table 1 and Fig- 
ure 1. The median spectral type of our sample is KG, the me- 
dian H magnitude is 5.75, the median distance is 22 pc, the 
median proper motion amplitude is 240 mas yr"', and the me- 
dian [Fe/H] is 0.00 dex (standard deviation of 0.21 dex). 

Despite our effort to select only single stars, our observa- 
tions show that 16 of the 85 target stars are close double or 
triple systems; this is indicated in the last column of Table 1. 
A thorough review of the literature revealed that 1 1 of these 
were known at the time the target list was compiled, two of 
which are astrometric multiples that had never been resolved 
prior to our observations (HD 14802 and HD 166181). Five 
other multiple systems were resolved with AO only after the 
target list was compiled (HD 77407, HD 129333, HD 135363, 
HD 160934, and HD 220140). Finally, the star HD 213845 is 
reported to be part of a binary system for the first time here. 
The multiple systems observed are discussed further in §4.3. 

Age estimates for the stars in our sample, needed to convert 
the observed contrasts into mass detection limits using evo- 
lution models of giant planets,^ are reported in Table 1 along 
with the references used for their determination. Whenever 
possible, we have used ages stated explicitly in the literature 
or the age of the association to which a star belongs. When 
no specific age estimate was available for stars taken from 
Wichmann et al. (2003), ages of 10-50 Myr or 50-150 Myr 
were assigned to the stars having a Uthium abundance above 
or comparable to that of the Pleiades, respectively. For other 
stars that have lithium and/or X-ray measurements, ages were 
estimated from a comparison of the Li I 6708 A equivalent 
width and/or the ratio of the X-ray to bolometric luminosity 
with Figures 3 and/or 4 of Zuckerman & Song (2004b) respec- 
tively. When lithium or X-ray measurements were not avail- 
able, the kinematic ages were used as lower limits while the 
ages derived from the chromospheric activity index, log/?HK, 
were used as upper limits, as Song et al. (2004) showed that 
the latter ages tend to be systematically higher than those de- 
rived from lithium abundance or X-ray emission. When only 
the value of log/?HK was available, the calibration of Donahue 
(1993)^ was used to obtain an age estimate. Finally, when 
only kinematics measurements were available for a given star, 
an age of 100-5000 Myr or 50-5000 Myr was assigned if 
the star is a possible member of the Local Association or the 
IC 2391 supercluster respectively. 

3. OBSERVATIONS AND IMAGE PROCESSING 

3.1. Data acquisition and observing strategy 

All observations were obtained at the Gemini North tele- 
scope with the Altair adaptive optics system (Herriot et al. 
2000) and the NlRl camera (Hodapp et al. 2003) (pro- 
grams GN-2004B-Q-14, GN-2005A-Q-16, GN-2005B-Q-4, 
GN-2006A-Q-5, and GN-2006B-Q-5). The //32 camera was 
used, yielding 0.022" pixel"' and a field of view of 22" x 22". 
The field lens of Altair, which improves the off-axis adap- 

^ It is assumed that any planet and its primary star would be coeval. 
' This cahbration is given expUcitly in Henry et al. (1996). 



live optics correction, was not used for any observation as it 
was not available for the first epoch observations. Because 
it introduces an undetermined field distortion, having used 
the field lens for the second epoch observations only would 
have complicated or prevented verification of the physical as- 
sociation of companion candidates identified in the first epoch 
observations. The observations were obtained in the narrow 
band filter CH4-short (1.54-1.65 /zm), for the following rea- 
son. According to evolution models (e.g. Baraffe et al. 2003), 
planetary mass objects older than 10-20 Myr should have an 
effective temperature below 1000 K. Because of the large 
amounts of methane and the increased collision induced ab- 
sorption by H2 in their atmosphere, the near-infrared /T-band 
flux of such objects is largely suppressed. It is thus more 
efficient to search for giant planets in either the / or the H 
band; the latter was preferred in this study because higher 
Strehl ratios are achieved at longer wavelengths. As the bulk 
of the 7/-band flux of cool giant planets is emitted in a nar- 
row band centered at ~1.58 /^m because of important absorp- 
tion by methane beyond 1.6 /^m, it is even more efflcient to 
search for these planets using the CH4-short filter, which is 
well matched to the peak of the emission. Based on evolution 
models and synthetic spectra of giant planets (Baraffe et al. 
2003), it is expected that the mean flux density of a planet in 
the NIRI CH4-short filter be between 1.5 and 2.5 times higher 
(0.44-1 .0 mag brighter) than in the broad band H filter, de- 
pending on the specific age and mass of the planet. These 
factors are consistent with the factors 1.6-2.0 (0.5-0.75 mag) 
calculated from the observed spectra of T7-T8 brown dwarfs, 
which have Jeff 800 K. 

The angular differential imaging (ADl, Marois et al. 2006) 
technique was used to suppress the PSF speckle noise and im- 
prove our sensitivity to faint companions. This technique con- 
sists of acquiring a sequence of many exposures of the target 
using an altitude/azimuth telescope with the instrument rota- 
tor turned off (at the Cassegrain focus) to keep the instrument 
and telescope optics aligned. This is a very stable configu- 
ration and ensures a high correlation of the sequence of PSF 
images. This setup also causes a rotation of the field of view 
(FOV) during the sequence. For each target image in such a 
sequence, it is possible to build a reference image from other 
target images in which any companion would be sufficiently 
displaced due to FOV rotation. After subtraction of the ref- 
erence image, the residual images are rotated to align their 
FOV and co-added. Because of the rotation, the residual PSF 
speckle noise is averaged incoherently, ensuring an ever im- 
proving detection limit with increasing exposure time. It has 
been shown that, for ADl with Altair/NIRI, the subtraction of 
an optimized reference PSF image from a target image can 
suppress the PSF speckle noise by a factor of ~12, and that 
a noise suppression factor of ^100 is achieved for the com- 
bination of 90 such difference images (Lafreniere et al. 2007; 
Marois et al. 2006). 

An individual exposure time of 30 seconds was chosen 
for all targets. This exposure time is long enough so that, at 
large separation, faint companion detection is limited by sky 
background noise rather than read noise, and short enough 
so that the radius of the region affected by saturation and 
non-linearity of the detector typically does not exceed 0.5". 
The nominal observing sequence consisted of 90 images, 
but oftentimes a few images had to be discarded due to brief 
periods of very bad seeing, loss of tracking, or the advent of 
clouds. No dithering was made during the main observing 
sequence to ensure a high correlation of the PSF images; flat- 
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TABLE 2 
GDPS OBSERVATION LOG 



Name 


Date 


Number of 


— - — — — 

Strehl 


— r : 

FOV rotation 


— ; 

Saturation 






exposures 


(%) 


(deg) 


j; „ ///\a 

radius ( ) 


HD 166 


2005/08/25 


83 


5-8 


55 


0.98 




2006/07/18 


83 


7-10 


81 


0.78 


HD 691 


2005/08/10 


90 


13-17 


70 


0.43 




2006/09/18 


117 


16-30 


88 


0.44 


HD 1405 


2004/08/22 


90 


4-10 


17 


0.53 




2005/08/04 


90 


6-18 


69 


0.40 


HD5996 


2005/08/12 


90 


18-20 


24 


0.50 




2006/09/25 


90 


15-17 


21 


0.50 


HD9540 


2005/08/14 


90 


16-19 


25 


0.55 




2006/09/28 


45 


14-17 


11 


0.61 


HD 10008 


2005/08/10 


90 


18-20 


36 


0.51 


GJ82 


2005/08/31 


90 


10-12 


27 


0.28 


HD 14802 


2005/08/20 


90 




23 


1.09 


HD 16765 


2005/09/10 


90 


14-17 


45 


0.72 


HD 17190 


2005/08/24 


90 


13-30 


108 


0.52 


HD 17382 


2004/12/22 


66 


15 


68 


0.55 




2005/09/1 1 


90 


19-23 


104 


0.52 


HD 17925 


2004/1 1/04 


83 


>18'' 


29 


0.66 


HD 18803 


2004/12/24 


90 


7-14 


99 


0.70 




2005/09/12 


78 


17-18 


108 


0.64 


HD 19994 


2005/08/3 1 


90 




44 


0.83 




2006/10/01 


57 




27 


0.72 


HD 20367 


2005/10/02 


90 


12-14 


67 


0.70 


2E 759 


2005/10/17 


59 


7-10 


31 


0.22 


HD 22049 


2005/09/08 


90 




32 


2.06 


HIP 17695 


2005/09/13 


89 


20-20 


45 


0.24 


HD 25457 


2005/10/02 


90 




43 


0.96 


HD 283750 


2004/10/24 


90 


15 


99 


0.54 




2005/10/04 


87 


19-23 


101 


0.59 


HD 30652 


2005/09/12 


52 


_ 


35 


1.86 


GJ182 


2004/11/05 


90 


16-20 


31 


0.37 




2005/10/17 


33 


11-11 


29 


0.39 


GJ 234A 


2005/11/05 


72 


16 


34 


0.42 


GJ281 


2005/03/25 


67 


9-10 


49 


0.52 




2006/02/12 


25 


8-9 


11 


0.47 


GJ285 


2005/03/18 


20 


_ 


10 


0.45 




2006/02/12 


90 


4-5 


73 


0.55 


HD 72905 


2005/04/23 


84 


7 


25 


0.87 


HD 75332 


2005/04/24 


89 


>17'' 


27 


0.59 




2006/12/20 


16 


>16b 


11 


0.59 


HD 77407 


2005/04/26 


84 


r6-19 


33 


0.61 


HD 78141 


2004/12/21 


85 


14-16 


19 


0.55 


HD 82558 


2005/04/18 


90 




30 


0.61 


HD 82443 


2004/12/25 


75 


18 


28 


0.61 


GJ 393 


2005/04/20 


90 


13-15 


44 


0.55 


HD 90905 


2005/03/18 


90 


13-18 


47 


0.61 




2006/04/11 


35 


13-15 


14 


0.55 


HD 91901 


2005/04/29 


71 


9 


22 


0.44 


HD 92945 


2005/05/26 


85 


15-16 


19 


0.61 




2006/05/16 


10 


10-11 


2 


0.53 


HD 93528 


2005/04/30 


86 




26 


0.39 


GJ402 


2005/04/26 


79 


12-16 


37 


0.39 




2006/02/16 


60 


6-10 


33 


0.35 


HD 96064 


2005/04/19 


89 


21-23 


37 


0.50 




2006/03/05 


90 


13-19 


36 


0.50 


HD 97334 


2005/04/18 


90 


16-17 


54 


0.70 


HD 102195 


2005/04/24 


91 


20-21 


54 


0.41 




2006/03/18 


82 


12-18 


30 


0.39 


HD 102392 


2005/04/23 


89 


19-24 


32 


0.39 




2006/03/12 


90 


9-13 


31 


0.40 


HD 105631 


2005/05/29 


90 


14-19 


45 


0.55 


HD 107146 


2005/05/30 


90 


21-26 


71 


0.57 


HD 108767B 


2005/04/22 


90 


14 


27 


0.43 




2006/02/16 


43 


10-11 


14 


0.41 


HD 109085 


2005/05/26 


90 




22 


1.09 




2006/03/12 


15 




3 


1.09 


BD+60 1417 


2005/04/18 


90 


18-23 


24 


0.26 




2006/04/1 1 


63 


12 


19 


0.24 


HD 111395 


2005/04/19 


89 


>12'' 


120 


0.77 


HD 113449 


2005/06/01 


47 


10-20 


37 


0.52 


GJ 507.1 


2005/06/07 


87 


5-7 


61 


0.44 


HD 116956 


2005/05/29 


90 


5-14 


27 


0.55 
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TABLE 2 — Continued 



Name 


Date 


Number of 


Strehl FOV rotation Saturation 






exposures 


(%) (deg) radius (")" 





2006/05/16 


60 


5-8 


18 


0.61 


HD 118100 


2005/04/27 


53 


_ 


18 


0.39 


GJ 524.1 


2005/04/18 


90 


18-25 


37 


0.26 




2006/05/18 


90 


13-14 


37 


0.22 


HD 124106 


2005/04/19 


86 


18-19 


32 


0.50 




2006/02/16 


80 


10-12 


24 


0.50 


HD 125161B 


2005/05/30 


90 


17-23 


31 


0.39 


HD 129333 


2005/04/20 


90 


19-20 


22 


0.48 


HD 130004 


2005/05/25 


90 


17-18 


105 


0.55 


HD 130322 


2005/05/27 


88 


15-19 


40 


0.48 




2006/05/15 


10 


10-10 


5 


0.40 


HD 130948 


2005/04/17 


90 


>9'' 


122 


0.83 


HD 135363 


2005/04/18 


87 


14-15 


19 


0.48 




2006/02/16 


60 


8-9 


14 


0.44 


HD 139813 


2005/05/30 


90 


>12'' 


20 


0.57 


HD 141272 


2005/04/19 


90 


r8-19 


47 


0.55 




2006/03/12 


42 


13 


20 


0.56 


HD 147379B 


2005/04/18 


90 


17-17 


22 


0.50 


GJ628 


2005/04/17 


90 


11 


29 


0.70 




2006/04/11 


40 


9-14 


13 


0.66 


HIP 81084 


2005/04/19 


73 


17-18 


30 


0.33 




2006/05/15 


90 


8-13 


31 


0.22 


HD 160934 


2005/04/18 


84 


17-24 


24 


0.35 




2006/09/17 


14 


12-14 


4 


0.34 


HD 162283 


2005/04/20 


120 


15-19 


45 


0.38 




2006/09/16 


100 


27-29 


31 


0.33 


HD 166181 


2005/04/17 


90 


16 


76 


0.59 




2006/09/18 


45 


18-21 


37 


0.48 


HD 167605 


2005/05/27 


90 


20 


22 


0.39 


HD 187748 


2005/05/25 


97 


15-19 


30 


0.66 




2006/09/15 


75 


>22^ 


21 


0.50 


GJ 791.3 


2005/05/26 


87 


9-19 


54 


0.42 


HD 197481 


2005/07/29 


68 


6-10 


21 


0.87 


HD 201651 


2005/06/27 


90 


18-23 


21 


0.38 




2006/09/14 


30 


19-21 


7 


0.38 


HD 202575 


2005/07/16 


90 


17-23 


75 


0.57 




2006/09/14 


30 


16-18 


9 


0.56 


GJ4199 


2004/08/23 


65 


10-13 


118 


0.44 




2005/08/04 


90 


15-23 


136 


0.39 


HD 206860 


2005/08/10 


34 


>13'' 


56 


0.77 




2006/06/26 


60 


>15'' 


80 


0.61 


HD 208313 


2005/06/27 


90 


23-23 


67 


0.46 




2006/06/25 


89 


14-22 


66 


0.55 


V383 Lac 


2005/07/26 


66 


13-17 


28 


0.42 




2006/06/30 


77 


15-18 


27 


0.32 


HD 213845 


2005/08/24 


90 




26 


0.81 




2006/07/06 


90 




24 


0.83 


GJ 875.1 


2005/08/10 


90 


16-18 


69 


0.33 




2006/07/07 


79 


7-17 


61 


0.31 


GJ876 


2005/08/21 


82 


9-16 


28 


0.68 


GJ 9809 


2005/08/04 


90 


18-20 


25 


0.31 




2006/09/14 


120 


25-27 


31 


0.22 


HD 220140 


2005/08/05 


90 


16-18 


21 


0.59 




2006/07/16 


82 


7-9 


19 


0.63 


HD 221503 


2005/08/31 


90 


21-22 


28 


0.52 


GJ900 


2004/08/24 


90 


15-21 


17 


0.46 




2005/09/08 


90 


16-22 


46 


0.42 


GJ 907.1 


2005/09/07 


65 


5-15 


22 


0.37 




2006/07/17 


44 


8 


16 


0.31 



^ Radius at which the PSF radial intensity profile reaches 75% of the detector well capacity. 

^ Only a lower estimate of the Strehl ratio can be obtained as the PSF peak is in the non-linear regime or sligthly saturated. 



field errors, bad pixels, and cosmic ray hits are naturally av- 
eraged/removed with ADI because of the FOV rotation. The 
PSF centroid was found to wander over the detector by typ- 
ically 2-5 pixels throughout an observing sequence because 
of mechanical flexure and differential refraction between the 
wavefront sensing and science wavelengths; for a handful of 
targets the variation sUghtly exceeded 10 pixels. Short unsat- 
urated exposures were acquired before and after the main se- 
quence of (saturated) images for photometric calibration and 
Strehl ratio estimation; these observations were acquired in 



sub-array mode (256 x 256 or 512x512 pixels), for which the 
minimum exposure time is shorter. Typically, an unsaturated 
sequence consisted of five exposures each obtained at a differ- 
ent dither position. The unsaturated observations are missing 
for a few targets as they were either skipped in the execution 
of the program, or they turned out to be saturated despite us- 
ing the shortest possible exposure time. Table 2 summarizes 
all observations. The last column of the table ("saturation ra- 
dius") indicates the separation at which the radial profile of 
the PSF reaches 75% of the detector full well capacity; linear- 
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Fig. 2. — Illustration of the ADI noise attenuation process. Panel (a) shows an original 30-s image of the young star HD 691 after subtraction of an azimuthally 
symmetric median intensity profile, panels (b) and (c) both show, with a different intensity scale, the corresponding residual image after ADI subtraction using 
the LOCI algorithm, and panel (d) shows the median combination of 117 such residual images. Display intensity ranges are ±5 X 10"* and ±10"* of stellar 
PSF peak for the top and bottom rows respectively. Each panel is 10" on a side. The diffraction spikes from the secondary mirror support vanes and the central 
saturated region are masked. The faint point source (Am = 14.9) visible in panel (d) at a separation of 2.43" and P.A. of 7.3° could not have been detected 
without ADI processing. 



ity should be better than 1% at this level (Hodapp et al. 2003). 
We have not analyzed the data inside this separation; point 
sources located at least one PSF fuU-width-at-half-maximum 
(FWHM) past this separation can be detected in our analysis, 
provided that their brigthness is above the detection limit. 

3.2. Data reduction 

For each sequence of short unsaturated exposures, a sky 
frame was constructed by taking the median of the images 
obtained at different dither positions; this sky frame was sub- 
tracted from each image. The images were then divided by 
a flat field image. The PSFs of a given unsaturated sequence 
were registered to a common center and the median of the im- 
age sequence was obtained. The center of the PSFs were de- 
termined by fitting a 2-dimensional Gaussian function. As an 
indication of the quaUty of an observing sequence, the Strehl 
ratio was calculated by comparing the peak pixel value of the 
observed PSF image with that of an appropriate theoretical 
PSF. The calculated Strehl ratio values are reported in Table 2; 



two values are indicated for a target when unsaturated data 
were obtained before and after the main saturated sequence. 
Strehl ratios were typically in the range 10-20%. 

Images of the main saturated sequence were first divided 
by a flat field image. Bad and hot pixels, as determined 
from analysis of the flat field image and dark frame respec- 
tively, were replaced by the median value of neighboring pix- 
els. Field distortion was corrected using an IDL procedure 
provided by the Gemini staff (C. Trujillo, private communi- 
cation) and modified to use the IDL interpolate function with 
cubic interpolation. The plate scale and field of view orien- 
tation for each image were obtained from the FITS header 
keywords. 

For each sequence of saturated images, the stellar PSF of 
the first image was registered to the image center by maxi- 
mizing the cross-correlation of the PSF diffraction spikes with 
themselves in a 180-degree rotation of the image about its 
center. The stellar PSF of the subsequent images was regis- 
tered to the image center by maximizing the cross-correlation 



10 



Lafreniere et al. 



of the PSF diffraction spikes with those in the first image. 
Prior to shifting, the 1024 x 1024 pixel images were padded 
with zeros to 1450 x 1450 pixel to ensure that no FOV would 
be lost. An azimuthally symmetric intensity profile was fi- 
nally subtracted from each image to remove the smooth see- 
ing halo. 

Next, the stellar PSF speckles were removed from each 
image by subtracting an optimized reference PSF image ob- 
tained using the "locally optimized combination of images" 
(LOCI) algorithm detailed in Lafreniere et al. (2007). The 
heart of this algorithm consists in dividing the target image 
into subsections and obtaining, for each subsection indepen- 
dently, an optimized reference PSF image consisting of a lin- 
ear combination of the other images of the sequence for which 
the rotation of the FOV would have displaced sufficiently an 
eventual companion. For each subsection, the coefficients of 
the linear combination are optimized such that its subtrac- 
tion from the target image minimizes the noise. The subsec- 
tions geometry and the algorithm parameters determined in 
Lafreniere et al. (2007) were used for all targets. The residual 
images were then rotated to align their FOV and their median 
was obtained. Figure 2 illustrates the PSF subtraction process. 

3.3. Photometric calibration and uncertainty 

As the stellar PSF peak is saturated for the main sequence of 
images, and since much image processing is done to subtract 
the stellar PSF from each image, special care must be taken 
to calibrate the photometry of the residual images and ensure 
that the contrast limits calculated are accurate. 

When the PSF peak is saturated, relative photometry can 
be calibrated by scaling the stellar flux measured in the un- 
saturated images obtained before and/or after the saturated 
sequence according to the ratio of the exposure times of the 
saturated and unsaturated images. However, the accuracy of 
this calibration method is affected by the (unknown) varia- 
tions in Strehl ratio, hence of the peak PSF flux, that may have 
occurred between the saturated and unsaturated observations. 
To mitigate this problem, the calibration approach we adopted 
relies on a sharp ghost artifact located (+0.09", -2.45") from 
the PSF center in the ALTAIR/NIRI images. Since the inten- 
sity of this ghost artifact is proportional to the PSF intensity, it 
can be used to infer the peak flux of a saturated PSF. This was 
verified for all sequences for which both unsaturated and satu- 
rated data were available. First, the stellar flux was measured 
in the unsaturated images using a circular aperture of diam- 
eter equal to the FWHM of the PSF. When unsaturated data 
were acquired both before and after the saturated sequence, 
the mean of the two values was used. Then the flux of the 
ghost artifact in the same aperture was measured for each im- 
age of the saturated sequence. The median of these values, 
scaled according to the ratio of the exposure times of the satu- 
rated and unsaturated images, was then compared to the stel- 
lar flux, and the process was repeated for all sequences that 
include both saturated and unsaturated data. Similar values 
were found for all sequences; the mean ratio of the flux of the 
ghost over that of the PSF peak was found to be 6.1 x 10"^, 
with a standard deviation of 0.6 x 10"^. Comparisons of the 
flux of background stars bright enough to be visible in each 
individual image of a sequence with the flux of the ghost in 
the corresponding images also confirmed that the intensity of 
the ghost is indeed directly proportional to the intensity of 
off-axis sources. 

The procedure used for calibrating the photometry was the 
following. The flux of the ghost was measured for each im- 
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Angular separation (arcsec) 

Fig. 3. — Typical values of /sub (solid line), /aniso (dotted line), and /sm 
(dashed line) as a function of angular separation. The curves shown are for 
the target HD 166181. 

age of a sequence and the median of these values, divided by 
the ratio quoted above, was taken to represent the peak stellar 
PSF flux, /v- This calibration method should be more accu- 
rate than the one based solely on unsaturated data obtained 
before and/or after the saturated sequence because the me- 
dian ghost flux is affected in the same way as the median of 
all residual images by the variations of Strehl ratio that may 
have occurred during the sequence of saturated images or be- 
tween the saturated and unsaturated measurements. For this 
reason, this calibration was used even for the sequences for 
which unsaturated data were available. 

Observations obtained with ALTAIR without the field lens 
suffer from important off-axis Strehl degradation because of 
anisoplanatism; this degradation must be taken into account 
when calculating contrast. Unfortunately, it is virtually im- 
possible to quantify the specific degradation pertaining to our 
data as there are no bright reference off-axis point sources 
available for every sequence of images. Instead, we have 
used the average anisoplanetism Strehl ratio degradation for- 
mula indicated on the ALTAIR webpage"*, which is /aniso(^) = 
5(6')/5o = where 5(9) is the Strehl ratio at angular 

separation 0, expressed in arcseconds, and So is the on-axis 
Strehl ratio. This factor was used to correct the noise and the 
flux of faint point sources measured in the residual images. 

As explained in Lafreniere et al. (2007), while the subtrac- 
tion of an optimized reference PSF obtained using the LOCI 
algorithm leads to better signal-to-noise (S/N) ratios, it re- 
moves partially the flux of the point sources sought after. This 
flux loss must be accounted for when calculating contrast. 
This is done by calculating the normalized residual intensity, 
/sub, of artificially implanted point sources after execution of 
the subtraction algorithm; the method used is described in 
§4.3 of Lafreniere et al. (2007). Then using flux measure- 
ments made in the residual image, the factor /sub is used to 
infer the true flux of a point source, i.e. that before execution 
of the subtraction algorithm. 

Another effect that must be taken into account for ADl 
data is the azimuthal smearing of an off-axis point source 
that occurs as the field of view rotates during an integration; 
this causes a fraction of the source's flux to fall outside of 
the circular aperture used for photometric measurements. 
The amount of flux loss in the aperture was calculated for each 

^ http://www.gemini.edu/sciops/instruments/altair/ 
altairCommissioningPerformance.html 
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TABLE 3 
Photometric uncertainties 



Sep. (") 


<4 


4-7 7-10 


10-13 > 13 


<T (mag) 


0.07 


0.12 0.15 


0.26 0.39 



sequence of images as follows. For a given angular separa- 
tion and for each image of a sequence, a copy of the unsatu- 
rated PSF was smeared according to its displacement during 
an integration. When unsaturated data were unavailable, a 2D 
Gaussian of the appropriate FWHM was used in place of the 
unsaturated PSF. The median of these smeared PSFs was ob- 
tained and the flux in a circular aperture was measured. This 
flux was divided by the flux of the original PSF in the same 
aperture to obtain the smearing factor f^^, which is used to 
correct the flux or noise measured in the images. 

Given all of these considerations, the contrast at angular 
separation 9 was calculated as 



/aniso(^)/sm(^)/sub(^) ^* ' 

where F(0) is either the noise or the flux of a point source 
in a circular aperture of diameter equal to one PSF FWHM, 
at angular separation 9, in the residual image. Note that the 
contrast in the equation above is defined such that a fainter 
companion, or a smaller residual noise, has a smaller con- 
trast value. Eq. (1) was used for all contrast calculations in 
the present work. Typical correction factors as a function of 
angular separation are shown in Figure 3. 

An estimate of the photometric accuracy resulting from the 
entire process was obtained by calculating the mean absolute 
difference between the magnitudes calculated at two epochs 
for every faint background star that was observed twice (see 
§4.2); this mean absolute difference was taken to represent 
■\/2 times the photometric uncertainty. This photometric un- 
certainty was found to vary significantly with angular sepa- 
ration, indicating that it is dominated by the uncertainty on 
the anisoplanatism factor. The photometric uncertainty as a 
function of angular separation is reported in Table 3; it is typ- 
icaUy 0.07-0.15 mag for separations below 10". For com- 
pleteness, it is noted that a higher photometric uncertainty, by 
about 0.08 mag, results when the unsaturated data obtained 
before and/or after the main sequence of saturated images are 
used to determine F^, rather than the median flux of the ghost 
artifact, justifying our choice to use the cahbration based on 
the flux of the ghost for all sequences. 

4. RESULTS 

4.1. Detection limits 

Detection limits are based on a measure of the noise in the 
residual images. To calculate this noise, the residual images 
were first convolved by a circular aperture of diameter equal 
to one PSF FWHM, which is typically ~0.07", and flie noise 
as a function of angular separation from the image center, 
F{9), was determined as the standard deviation of the pixel 
values in an annulus of width equal to one PSF FWHM. As 
shown in Lafreniere et al. (2007) and Marois et al. (2007b), 
the noise in an ADI residual image has a distribution similar 
to a Gaussian; using a 5cr detection threshold is thus appropri- 
ate for our data to limit the number of false positives. Given 
that a residual image typically contains ^2x10^ resolution 



elements, roughly 0. 1 false positive per target is expected on 
average. Because of the underlying noise in a residual im- 
age, some sources near the detection threshold might not be 
detected. From Gaussian statistics, the probability that the 
residual signal underlying a source is below Oct, -Ict, or -2ct 
is 50%, 16%, or 2.3%, respectively. Our detection complete- 
ness for sources whose true intensities are 5ct, 6a, or >7ct is 
thus 50%, 84%, or >97.7%, respectively. Note that for a sim- 
ilar reason, some sources whose true intensities are below the 
5ct threshold could be detected as well. These effects wiU be 
taken into account appropriately in the statistical analysis of 
the results presented in §5. 

The detection limits achieved for all target stars, expressed 
in magnitude difference, are presented in Table 4. The last 
two Unes of this table present the median and best contrast, 
over the 85 observations, achieved at each angular separation. 
The median detection limits in magnitude difference are 9.5 
at 0.5", 12.9 at 1", 15.0 at 2", and 16.5 at 5". The detec- 
tion limits are presented graphically in Figure 4 for the stars 
HD 208313, HD 166181, andGJ 507.1, which are representa- 
tive of poor, median, and good contrast performance, respec- 
tively. 

For consistency we have verified the validity of these detec- 
tion limits by implanting fiducial sources in the sequence of 
original images and then processing the data as described in 
§3.2. An example, incorporating artificial sources at the 5ct 
and IOct levels at various separations, is shown in Figure 5 for 
the stars HD 208313, HD 166181, and GJ 507.1. As visible in 
this figure, sources exactly at our detection limits can indeed 
be detected with the expected completeness level. 

One must resort to evolution models of giant planets to 
convert the detection limits mentioned above into mass 
limits. Traditionally, such evolution models have assumed 
arbitrary initial conditions for the planets (e.g. Baraffe et al. 
2003; Burrows et al. 1997), with the caution that their results 
depend on the specific initial conditions adopted for ages 
below a few million years (Baraffe et al. 2002). Recent 
evolution models (Marley et al. 2007) that incorporate 
initial conditions calculated exphcitly for planets formed 
through core accretion indicate that it may in fact take as 
much as 10-100 Myr before the planets "forget" their initial 
conditions; the effect being more important for more massive 
planets. Nevertheless, given the typical ages of our target 
stars (50-300 Myr) and the good contrast limits we have 
reached, the different evolution models should yield similar 
mass detection limit estimates. As a simple example, consider 
a contrast of 12.9 mag in the NlRl CH4-short filter around a 
KO star (typical at a separation of 1"). The "hot start" models 
of Baraffe et al. (2003) would give masses of 2.6 Mjup and 
3.9 Mjup at 50 Myr and 100 Myr, respectively, while the "core 
accretion" models of Marley et al. (2007) would give masses 
of ^3.0 Mjup and ^4.5 Mjup, respectively^. The difference 
between the models would be smaller for smaller masses 
(better contrast limits (i.e. beyond ^1") and/or greater ages), 
while it would be larger for larger masses (worse contrast 
limits and/or smaUer ages). In this work, keeping the latter 
caveat in mind, we have used the COND evolution models of 
Baraffe et al. (2003), for which absolute H-band magnitudes 
as a function of mass and age are readily available. The 
following procedure was used to estimate the contrast, in the 



' For this simple calculation, it was assumed that the luminosity ratios 
between the "hot start" and "core accretion" models were representative of 



the //-band magnitude differences. 



12 



Lafreniere et al. 



Projected physical separation (AU) at 22 pc 



25. 50. 75. 100. 125. 150. 175. 200. 
8| ' I X 




181 , , , I , , , I , , , I , , , I , , , I 

2 4 6 8 10 

Angular separation (arcsec) 



Fig. 4. — Survey detection limits in difference of magnitude (in the NIRI CH4-short filter) between an off-axis point source and the target star, at the 6cr level. 
The top, middle, and bottom curves are respectively for the targets GJ 507.1, HD 166181, and HD 208313, which are representative of poor, median, and good 
performance reached by the survey. Companion candidates identified around targets of galactic latitude \b\ < 15 are shown by + symbols, while those identified 
around targets with \b\ > 15 are shown by X symbols. The two filled circles near (2.6,8.6) indicate the components of the binary brown dwarf companion to 
HD 130948. The fiducial point sources shown in Fig. 5 are marked with triangles. The top and right axes show, for reference only, the projected separation in AU 
and the detection limits in Mj^p that would apply for a 100 Myr old KO star located 22 pc away. 



TABLE 4 

GDPS DETECTION LIMITS' 



Name 


0.50" 


0.60" 


0.75" 


1.00" 


1.25" 


1.50" 


2.00" 


2.50" 


3.00" 


4.00" 


5.00" 


7.50" 


10.00" 


HD 166 








12.5 


13.1 


13.9 


14.9 


15.4 


15.9 


16.5 


16.9 


17.3 


17.3 


HD691 




11.1 


12.1 


13.2 


14.1 


14.7 


15.6 


15.9 


16.2 


16.6 


16.7 


16.6 


16.3 


HD 1405 


9.2 


10.5 


11.4 


12.7 


13.5 


14.0 


14.8 


15.3 


15.7 


16.0 


16.1 


16.1 


15.8 


HD 5996 




10.8 


12.0 


13.2 


14.1 


14.6 


15.4 


15.8 


16.1 


16.5 


16.6 


16.6 


16.4 


HD 9540 






11.8 


13.1 


14.0 


14.5 


15.4 


16.0 


16.4 


17.0 


17.3 


17.6 


17.5 


HD 10008 




10.0 


11.2 


12.4 


13.2 


13.8 


14.7 


15.2 


15.6 


16.2 


16.5 


16.5 


16.3 


GJ82 


8.9 


9.5 


10.5 


11.8 


12.5 


13.2 


13.7 


14.3 


14.6 


14.9 


15.0 


14.8 


14.6 


HD 14802 










11.8 


12.4 


13.3 


14.0 


14.7 


15.8 


16.8 


17.4 


17.9 


HD 16765 








13.0 


13.9 


14.5 


15.3 


15.8 


16.2 


16.9 


17.4 


17.5 


17.6 


HD 17190 




10.5 


12.2 


13.7 


14.2 


14.8 


15.5 


15.9 


16.3 


16.6 


16.8 


16.6 


16.2 


HD 17382 




10.8 


12.0 


13.3 


14.1 


14.6 


15.4 


15.9 


16.3 


16.8 


17.0 


17.0 


16.7 


HD 17925 






11.9 


13.6 


14.6 


15.4 


16.2 


16.8 


17.1 


17.6 


17.7 


17.7 


17.4 


HD 18803 






11.3 


12.9 


13.8 


14.5 


15.5 


16.0 


16.5 


16.8 


17.1 


17.2 


16.9 


HD 19994 








13.5 


14.3 


15.0 


15.8 


16.4 


16.7 


17.4 


17.8 


18.3 


18.4 


HD 20367 








11.6 


12.2 


12.8 


13.9 


14.4 


14.8 


15.6 


16.0 


16.3 


16.1 


2E 759 


8.6 


9.4 


9.9 


11.0 


11.8 


12.2 


13.0 


13.4 


13.6 


13.9 


14.0 


13.9 


13.6 


HD 22049 
















15.9 


16.5 


17.3 


17.7 


18.5 


18.9 


HIP 17695 


10.0 


10.8 


11.8 


12.8 


13.6 


14.2 


14.8 


15.1 


15.3 


15.7 


15.7 


15.6 


15.3 


HD 25457 










12.5 


13.1 


13.9 


14.8 


15.2 


16.0 


16.6 


17.0 


17.0 


HD 283750 






12.2 


13.4 


14.2 


15.1 


15.9 


16.4 


16.8 


17.2 


17.2 


17.1 


16.7 


HD 30652 














14.9 


15.5 


15.9 


16.7 


17.3 


18.2 


18.6 


GJ 182 


10.0 


10.5 


11.9 


13.1 


14.0 


14.7 


15.4 


15.8 


16.1 


16.4 


16.5 


16.4 


16.2 


GJ 234A 


9.5 


10.1 


11.2 


12.3 


13.3 


13.9 


14.6 


15.1 


15.4 


15.9 


16.2 


16.3 


16.1 


GJ281 




9.0 


10.4 


12.0 


12.9 


13.5 


14.3 


14.6 


15.0 


15.3 


15.3 


15.4 


15.2 


GJ285 




8.0 


10.1 


11.6 


12.6 


13.3 


13.8 


14.5 


14.9 


15.5 


15.8 


15.9 


15.8 


HD 72905 








11.2 


12.5 


13.1 


14.2 


14.9 


15.4 


16.3 


16.7 


17.4 


17.7 


HD 75332 






10.8 


12.3 


13.0 


13.9 


14.9 


15.5 


15.7 


16.6 


17.1 


17.4 


17.3 


HD 77407 






10.3 


11.4 


12.3 


13.0 


14.0 


14.8 


15.0 


15.7 


16.0 


16.3 


16.2 


HD 78141 






11.5 


13.0 


13.7 


14.5 


15.4 


15.8 


16.1 


16.5 


16.6 


16.5 


16.3 


HD 82558 






11.5 


12.9 


13.8 


14.4 


15.4 


15.9 


16.1 


16.6 


16.8 


17.0 


16.7 


HD 82443 






11.5 


13.0 


14.1 


14.8 


15.9 


16.4 


16.8 


17.2 


17.5 


17.7 


17.5 
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Name 


0.50" 


0.60" 


0.75" 


1 .00" 


1.25 


1.50 


2.00" 


2.50" 


3.00" 


4.00" 


5.00" 


7.50" 


10.00" 


GJ 393 






11.8 


13.3 


14.1 


14.6 


15.6 


16.0 


16.2 


16.7 


16.8 


16.9 


16.8 


HD 90905 






11.4 


12.7 


13.7 


14.1 


15.1 


15.7 


16.1 


16.5 


16.6 


16.6 


16.4 


HD 91901 




9.2 


10.0 


11.4 


12.1 


12.8 


13.6 


14.1 


14.4 


14.9 


14.8 


14.8 


14.6 


HD 92945 






10.8 


12.1 


13.0 


13.8 


14.6 


15.1 


15.5 


15.9 


16.1 


16.3 


16.1 


HD 93528 


8.5 


9.3 


10.2 


11.6 


12.6 


13.3 


14.2 


14.8 


15.0 


15.5 


15.7 


15.9 


15.7 


GJ402 


8.4 


9.2 


10.5 


11.6 


12.5 


13.1 


14.0 


14.5 


14.9 


15.4 


15.4 


15.6 


15.3 


HD 96064 




10.9 


12.3 


13.5 


14.3 


14.9 


15.6 


16.1 


16.3 


16.6 


16.8 


16.8 


16.6 


HD 97334 








13.6 


14.7 


15.1 


16.0 


16.4 


16.7 


17.2 


17.4 


17.5 


17.3 


HD 102195 


9.8 


11.2 


12.2 


13.3 


14.1 


14.7 


15.4 


15.9 


16.1 


16.5 


16.6 


16.6 


16.3 


HD 102392 


9.5 


10.3 


11.4 


12.6 


13.5 


13.9 


14.7 


15.3 


15.6 


16.1 


16.3 


16.3 


16.2 


HD 105631 






11.7 


12.8 


13.5 


14.2 


15.1 


15.5 


16.0 


16.4 


16.7 


16.8 


16.5 


HD 107146 






11.7 


12.5 


13.5 


14.0 


15.0 


15.4 


15.8 


16.2 


16.5 


16.5 


16.3 


HD 108767B 


8.4 


9.7 


10.6 


11.9 


12.8 


13.5 


14.3 


14.9 


15.1 


15.7 


15.8 


16.0 


15.7 


HD 109085 










13.4 


14.0 


14.9 


15.8 


16.3 


17.2 


17.7 


18.3 


18.5 


BD+60 1417 


10.0 


11.1 


12.0 


13.0 


13.8 


14.2 


14.7 


15.0 


15.3 


15.5 


15.5 


15.4 


15.1 


HD 111395 








13.4 


14.3 


15.0 


15.9 


16.4 


16.7 


17.2 


17.4 


17.6 


17.3 


HD 113449 






11.5 


12.6 


13.7 


13.9 


14.9 


15.4 


15.7 


16.3 


16.5 


16.6 


16.4 


GJ 507.1 




9.6 


10.4 


11.5 


12.2 


12.9 


13.7 


14.3 


14.6 


15.0 


15.2 


15.2 


14.9 


HD 1 16956 






11.3 


12.7 


13.5 


14.2 


15.1 


15.7 


16.0 


16.5 


16.7 


16.8 


16.6 


HD 118100 


8.4 


9.4 


10.5 


11.6 


12.3 


12.8 


13.5 


14.0 


14.1 


14.4 


14.5 


14.4 


14.2 


GJ 524.1 


10.1 


1 1.0 


12.0 


13.0 


13.6 


14.2 


14.9 


15.2 


15.4 


15.4 


15.5 


15.4 


15.0 


HD 124106 




10.3 


11.6 


13.0 


13.8 


14.4 


15.4 


15.7 


15.9 


16.4 


16.7 


16.8 


16.6 


HD 125161B 


10.5 


11.3 


12.4 


13.6 


14.3 


14.6 


15.4 


15.8 


16.0 


16.2 


16.4 


16.3 


16.1 


HD 129333 




10.7 


11.7 


13.2 


13.9 


14.4 


15.3 


15.7 


16.2 


16.4 


16.7 


16.7 


16.5 


HD 130004 






12.0 


13.1 


14.1 


14.5 


15.3 


15.8 


16.1 


16.5 


16.7 


16.7 


16.4 


HD 130322 




11.1 


12.1 


13.2 


13.9 


14.3 


15.2 


15.6 


15.9 


16.3 


16.4 


16.4 


16.2 


HD 130948 








12.4 


13.2 


13.8 


14.7 


15.4 


15.7 


16.5 


16.9 


17.3 


17.3 


HD 135363 




9.2 


10.9 


12.3 


13.1 


13.7 


14.6 


15.1 


15.3 


15.5 


15.7 


15.6 


15.4 


HD 139813 






10.3 


11.2 


11.9 


12.6 


13.6 


14.3 


14.9 


15.7 


16.1 


16.3 


16.1 


HD 141272 






12.2 


13.7 


14.4 


15.0 


15.8 


16.2 


16.5 


16.9 


16.9 


17.1 


16.9 


HD 147379B 




10.0 


11.3 


12.8 


13.5 


14.1 


15.0 


15.3 


15.6 


15.8 


16.0 


16.0 


15.7 


GJ 628 






10.4 


12.2 


13.0 


13.7 


14.6 


15.2 


15.7 


16.2 


16.6 


16.9 


16.7 


HIP 81084 


9.5 


10.3 


11.4 


12.3 


13.0 


13.5 


14.0 


14.4 


14.6 


14.7 


14.7 


14.6 


14.3 


HD 160934 


9.5 


10.1 


11.2 


12.5 


13.3 


13.9 


14.6 


14.9 


15.0 


15.3 


15.3 


15.2 


14.9 


HD 162283 


10.3 


11.2 


12.2 


13.4 


14.0 


14.6 


15.2 


15.7 


16.1 


16.4 


16.5 


16.5 


16.1 


HD 166181 




10.8 


11.7 


13.0 


13.7 


14.3 


15.0 


15.4 


15.8 


16.2 


16.5 


16.5 


16.3 


HD 167605 


9.4 


10.5 


11.4 


12.5 


13.3 


14.0 


14.8 


15.1 


15.6 


15.9 


16.1 


16.1 


15.9 


HD 187748 




10.8 


11.7 


12.9 


13.7 


14.5 


15.3 


15.9 


16.3 


17.0 


17.3 


17.6 


17.4 


GJ 791.3 


9.6 


11.0 


12.0 


13.3 


13.8 


14.4 


15.1 


15.6 


15.7 


16.0 


16.1 


16.1 


15.7 


HD 197481 








11.0 


11.7 


12.4 


13.5 


14.3 


14.7 


15.5 


16.1 


16.4 


16.3 


HD 201651 


10.1 


11.4 


12.3 


13.3 


14.1 


14.6 


15.3 


15.8 


16.1 


16.4 


16.5 


16.5 


16.3 


HD 202575 






11.4 


12.5 


13.3 


14.0 


14.9 


15.5 


16.1 


16.6 


16.8 


17.0 


16.7 


GJ 4199 


10.5 


11.2 


12.0 


13.2 


13.8 


14.5 


15.1 


15.6 


15.7 


16.0 


16.1 


15.8 


15.4 


HD 206860 






12.2 


13.3 


13.8 


14.5 


15.2 


15.7 


16.0 


16.5 


16.9 


17.0 


16.7 


HD 208313 




11.9 


13^0 


14.0 


UJ 


15^2 


16^0 


16^5 


16J 


17.2 


17^3 


17!3 


n.o 


V383 Lac 


10.2 


11.0 


11.9 


13.0 


13.6 


14.3 


14.8 


15.2 


15.5 


15.9 


16.1 


16.0 


15.8 


HD 213845 


- 


- 


- 


13.3 


14.0 


14.7 


15.7 


16.3 


16.8 


17.2 


17.6 


18.1 


18.0 


GJ 875.1 


9.6 


10.5 


11.2 


12.3 


13.1 


13.5 


14.4 


14.9 


15.1 


15.5 


15.6 


15.5 


15.1 


GJ876 








11.0 


12.2 


12.6 


13.7 


14.3 


15.1 


15.8 


16.2 


16.6 


16.6 


GJ 9809 


11.3 


12.1 


12.8 


14.0 


14.6 


15.0 


15.5 


15.9 


15.9 


16.2 


16.3 


16.1 


15.8 


HD 220140 






12.0 


13.1 


13.9 


14.5 


15.3 


15.8 


16.1 


16.4 


16.6 


16.5 


16.3 


HD 221503 




10.4 


11.8 


13.2 


14.1 


14.6 


15.3 


15.8 


16.2 


16.7 


17.0 


17.1 


17.0 


GJ900 


8.9 


10.1 


10.8 


12.4 


13.2 


13.9 


14.9 


15.4 


15.8 


16.1 


16.2 


16.1 


16.0 


GJ 907.1 


8.4 


9.0 


10.0 


11.2 


12.1 


12.5 


13.4 


14.0 


14.3 


14.8 


15.0 


15.1 


14.9 


Median 


9.5 


10.5 


11.5 


12.9 


13.6 


14.2 


15.0 


15.5 


15.8 


16.3 


16.5 


16.5 


16.3 


Best 


11.3 


12.1 


13.0 


14.0 


14.7 


15.4 


16.2 


16.8 


17.1 


17.6 


17.8 


18.5 


18.9 



^ Magnitude difference in the NIRI CH4-short filter, at a 5ct level. 

et al. 2003; AUard et al. 2001)^ We recall here (c.f. §3) 
that the ratio is typically 1.5-2.5 for giant planets de- 
pending on their mass and age. The stellar magnitudes in the 
NIRI CH4-short and broad band H filters were assumed to 
be equal, such that the contrast of the planet was obtained as 
nipi - Hi,, where //* is the //-band apparent magnitude of the 
target star. The 5 a contrast levels of planets of various masses 
orbiting a KO primary of 100 Myr,^ the typical target of the 
survey, are presented in Figure 4. For a typical target located 
at 22 pc from the Sun, the median detection Umits correspond 

* Spectra available at ftp://ftp.ens-lyon.fr/pub/users/CRAL/fallard/ 
^ An //-band absolute magnitude of 4.0 was used, tliis is the mean value 
of the KO stars in the sample. 



NIRI CH4-short filter, of a planet of given mass orbiting a 
given target. The absolute //-band magnitude of the planet 
was first obtained directly from the evolution models of 
Baraffe et al. (2003) and converted into an apparent magni- 
tude, //pi, using the distance of the star. The corresponding 
magnitude in the NIRI CH4-short filter was then calculated as 

mpi=//pi-2.51og(^^^, (2) 

where fcm and fn are the mean flux density of the planet in 
the NIRI CH4-short and broad band H filters, respectively; 
their values were calculated using a synthetic spectrum of ap- 
propriate effective temperature and surface gravity (Baraffe 
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Fig. 5. — Final S/N residual images for three sequences of images to which fiducial point sources have been implanted. The fiducial point sources have been 
added at 5 P.A.'s (0°, 72°, 144°, 216°, and 288°) and 3 angular separations (0.6", 1.0", and 2.0"). For bottom panels (d-f) the intensity of each source was set 
to the corresponding detection limit (5cr) indicated in Table 4, while it was set 0.75 mag brighter (i.e. lOcr) for top panels (a-c). Panels (a,d), (b,e), and (c,f) are 
for the stars HD 208313, HD 166181, and GJ 507.1, respectively. The bright spot at the upper left corner of panels (a,d) is a real background star. The display 
intensity scale is linear from -2 to +10 for top panels (a-c), and from —1 to +5 for bottom panels (d-f). In bottom panels (d-f), the sources that would have been 
detected (S/N>5) have been circled in white. According to expectations, the detection completeness is roughly 50% for sources whose true intensity is equal to 
the 5(T detection limit. 



to 10.8 Mjup at 1 1 AU, 3.9 Mj^p at 22 AU, 1.9 Mjup at 44 AU, 
and 1.4Mjup at 110 AU. 

The typical contrast reached by our survey improves on ear- 
lier surveys (e.g. Lowrance et al. 2005; Masciadri et al. 2005; 
Chauvin et al. 2006; Biller et al. 2007) by at least 1 mag at 
1", 1.5 mag at 2", and ^ 3 mag at larger separations. For 
the 27 targets for which our data were in the linear regime of 
the detector at a separation of 0.5", our detection hmits at this 
separation are similar to those achieved with the SDI device 
at the Very Large Telescope (Biller et al. 2007). The contrast 
reached by GDPS observations is the highest that has been 
achieved to date at separations larger than ~0.7". 

4.2. Candidate companion detections 

To identify candidate companions, the residual images were 
first convolved by a circular aperture of diameter equal to one 
PSF FWHM, and then converted to signal-to-noise (S/N) im- 
ages that were visually inspected for point sources at a >5(t 
level. After identification of a point source, its position was 
measured by fitting a 2D Gaussian function, and its flux was 
measured in an aperture of diameter equal to one PSF FWHM; 
both operations were done in the non-convolved residual im- 
age. The contrast of the point source was then calculated us- 
ing Eq.(l). More than 300 faint point sources were found 
around 54 targets, 188 of which are found around only 7 stars 
located at low galactic latitudes {b < 11°). Up to now, all but 
six of the 54 stars with candidates were re-observed at a sub- 
sequent epoch to verify whether or not the faint point sources 
detected are co-moving with the target star. 

AU candidate exoplanets observed at two epochs have been 



2.52 




fi r : 

0.0 0.2 0.4 0.6 0.8 1.0 

Epoch - 2005.6072 



Fig. 6. — Verification of the background nature of the point source detected 
aroimd the yoimg star HD 691. Open diamonds mark the observed separation 
(top) and P.A. {bottom) of the point source at the two epochs. The solid line 
indicates the expected separation and P.A. of a distant background source 
as a function of time. The observations agree very well with the expected 
motion of a background source, indicating that the source is not associated 
with HD 691. 

confirmed to be background sources by comparing their dis- 
placement between the two epochs with the expected dis- 
placement of a distant background source, based on the proper 
motion and parallax of the target; an example of this verifica- 
tion is presented in Figure 6. As a reference for future planet 
searches, a compilation of all faint point sources identified 
around our target stars is presented in Table 5. 
An estimate of the uncertainties on the measured separa- 
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tions and P. A. was obtained by calculating the mean abso- stars, the differences observed are dominated by our measure- 
lute difference between the separation and PA. measured at ment uncertainties. The mean absolute differences calculated 
the second epoch and those predicted for this epoch based on are taken to represent a/2 times the uncertainties; values of 
the parallax and proper motion of the target stars. Given the agep = 0.015" and ctp.a. = 0.2° are found, 
high precision on the parallax and proper motion of the target 



TABLE 5 
Point sources detected 



Star 


Epoch 


Separation^ 


FA.* 


Am'' 






(arcsec) 


(deg) 


mag 


HD 166 


2005.6482 


10.23 


82.9 


12.60 


HD 691 


2005.6072 


2.49 


12.1 


14.91 


HD 1405 


2004.6409 


3.95 


254.0 


13.98 


HD 5996 


2005.6128 


2.98 


118.6 


12.51 




2UU5.0125 


4.78 


71.6 


13.23 




onAC 1 OS 
ZUUj.OiZo 


J.OO 


OAS O 
ZOo.V 


1 ^ OQ 

1 j.vy 






6.95 


73.9 


15.43 




2005.6128 


9.11 


280.8 


15.72' 




2005.6128 


9.15 


228.1 


15.67 




2005.6128 


9.50 


120.6 


15.85' 




2005.6128 


9.58 


205.2 


13.51 




2005.6128 


9.94 


355.6 


13.86 




2005.6128 


10.33 


221.6 


10.86 




2005.6128 


10.49 


320.2 


10.93 




2005.6128 


10.55 


296.5 


14.41 




ZUUj.OlZiS 


11.10 


91 ft Q 






9005 61 ?S 


13.09 


215.7 


14.63 




2005.6128 


14.05 


141.1 


13.35 




2005.6128 


14.57 


314.9 


11.48 




2005.6128 


15.06 


137.5 


11.48 


HD9540 


2005.6182 


5.51 


308.8 


14.63 




2005.6182 


6.70 


120.9 


15.56 


GJ 82 


2005.6647 


4.24 


78.2 


13.89'' 






5 45 


35 3 


13 46'' 




zUUj.t)D4 / 


6.27 


157.0 


1 1 .UO 




2005.6647 


6.29 


228.1 


13.80'' 




2005.6647 


6.42 


307.0 


13.08'' 




2005.6647 


6.75 


105.7 


9.07'' 




2005.6647 


6.83 


106.4 


9.23'' 




2005.6647 


6.95 


25.7 


6.57'' 




2005.6647 


7.63 


117.2 


13.34'' 




2005.6647 


8.82 


315.9 


9.27'' 




2005.6647 


9.68 


318.7 


13.59'' 




2005.6647 


9.74 


334.0 


13.11'' 




2005.6647 


11.37 


228.3 


12.52'' 


HD 17382 


2004.9740 


11.78 


130.8 


13.16 


HD 18803 


2004.9795 


7.61 


166.1 


17.00 




2004.9795 


7.98 


208.3 


15.68 




2004.9795 


10.36 


52.8 


15.10 


HD 19994 


2005.6648 


6.18 


187.4 


17.62 




2005.6648 


6.30 


185.3 


16.06 




2005.6648 


11.64 


72.7 


17.51 


HD 283750 


2004.8132 


7.73 


175.8 


14.65 




onA/1 s 1 TO 


iZ. /Z 


1 r\A 
iu^.z 






ZUUj.OV / o 


Z.U4 


1 HA 1 


1 ^ 1 sd 
i J. 1 o 




ZUUJ.OV / o 


Q ^'i 


lA 1 A 
Z41 .4 


1 ft 


GJ 182 


9004. S45Q 


5.15 


220 3 


12 80 




2004.8459 


7.44 


233.7 


10.61 


GJ 234A 


2005.8455 


3.27 


48.8 


13.71'' 




2005.8455 


6.64 


304.9 


16.01'' 




2005.8455 


7.45 


215.1 


15.03'' 




2005.8455 


10.08 


179.3 


13.36'' 




2005.8455 


10.24 


84.3 


10.46'' 




2005.8455 


11.75 


103.0 


12.58'' 


GJ281 


2005.2286 


5.74 


237.0 


12.48 




2005.2286 


8.80 


288.4 


12.92 




2006.1158 


10.64 


224.6 


13.43f 


GJ285 


2005.2095 


8.83 


114.3 


11.55 


HD 75332 


2005.3107 


8.25 


141.7 


11.49 


HD 82443 


2004.9829 


5.27 


190.3 


11.64'' 




2004.9829 


5.42 


191.5 


16.65'' 




2004.9829 


8.33 


97.8 


12.24'' 




2004.9829 


10.17 


164.8 


le.w 




2004.9829 


13.74 


215.0 


14.67'' 
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Star 


Epoch 


Separation" 

(arcsec) 


RA.* 

(deg) 


Am'' 

mag 


HD 90905 


2005, 


.2098 


5.47 


188.2 


10.91 




2005, 


.2098 


12.41 


176.8 


13.32 


HD 92945 


2005 


.3983 


9.77 


236.2 


12.82 


HD 93528 


2005 


.3271 


4.82 


332.3 


14.27'' 


GJ 402 


2006, 


.1273 


12.46 


324.0 


10.80f 




2005 


.3164 


13.88 


337.7 


10.10 


HD 96064 


2005, 


.2972 


4.69 


29.7 


15.99"= 




2005, 


.2972 


5.94 


213.8 


8.18 




2005, 


.2972 


6.11 


213.6 


10.76 




2005, 


.2972 


8.90 


329.7 


15.90 


HD 102195 


2005, 


.3109 


11.94 


185.5 


13.82 


HD 102392 


2005 


.3082 


5.72 


42.0 


15.02 




2005 


.3082 


10.57 


308.9 


14.99= 


HD 108767B 


2005 


.3055 


6.72 


87.7 


12.41 




2005 


.3055 


8.28 


100.1 


14.62 




2005 


.3055 


10.20 


123.9 


15.10 


HD 109085 


2005 


.3986 


12.92 


256.2 


15.80 


BD+60 1417 


2005, 


.2946 


2.05 


298.4 


8.76 




2005 


.2946 


14.08 


133.5 


12.80"= 


HD 116956 


2005 


.4067 


9.34 


17.4 


15.05 


GJ 524.1 


2005, 


.2948 


7.59 


19.7 


13.09 


HD 124106 


2005, 


.2975 


7.51 


124.8 


13.85 




2005, 


.2975 


9.39 


342.2 


9.51 




2005, 


.2975 


9.60 


341.1 


8.71 




2005 


.2975 


10.39 


287.5 


14.65 




2005, 


.2975 


11.17 


291.7 


14.18 




2005, 


.2975 


12.06 


120.6 


15.45 


HD 130322 


2005, 


.4014 


7.61 


329.8 


11.00 


HD 130948 


2005 


.2922 


2.60 


103.1 


8.568 




2005 


.2922 


2.66 


104.0 


8.838 


HD 135363 


2005 


.2949 


7.50 


122.1 


10.28 


HD 139813 


2005 


.4097 


6.85 


271.3 


14.48'' 




2005 


.4097 


7.36 


272.2 


14.79'' 


HD 141272 


2005 


.2977 


2.31 


12.3 


11.44 




2005 


.2977 


4.03 


286.9 


16.59= 




2005 


.2977 


8.04 


305.4 


15.78= 




2005 


.2977 


8.32 


258.1 


15.66 




2005 


.2977 


10.95 


299.2 


9.90 




2005 


.2977 


11.43 


190.4 


15.14 




2005 


.2977 


12.26 


209.5 


12.31 


GJ628 


2005, 


.2924 


5.19 


259.0 


16.32"= 




2005, 


.2924 


6.29 


161.6 


14.81 




2005, 


.2924 


10.25 


232.1 


14.90 




2005, 


.2924 


10.52 


2.8 


13.42"= 




2005, 


.2924 


11.49 


308.6 


14.88 




2005, 


.2924 


12.72 


228.6 


15.10 




2005, 


.2924 


13.64 


215.4 


15.03 


HIP 81084 


2005, 


.2979 


6.84 


234.6 


13.10 




2005 


.2979 


8.64 


4.5 


13.70 




2005 


.2979 


9.69 


49.6 


11.22 




2005, 


.2979 


9.98 


68.4 


13.04 




2005 


.2979 


14.04 


226.4 


11.60 


HD 160934 


2005 


.2952 


4.08 


319.2 


11.03 




2005 


.2952 


8.94 


232.9 


9.96 


HD 162283 


2005 


.3007 


2.76 


118.7 


15.54 




2005 


.3007 


3.41 


154.4 


14.13 




2005 


.3007 


3.26 


4.9 


15.21 




2005, 


.3007 


3.68 


244.8 


14.01 




2005, 


.3007 


3.94 


152.8 


14.22 




2005, 


.3007 


4.22 


158.8 


12.52 




2005, 


.3007 


4.47 


111.9 


15.21 




2005, 


.3007 


4.37 


23.1 


13.23 




2005, 


.3007 


4.70 


299.8 


15.44 




2005, 


.3007 


4.95 


80.0 


10.82 




2005, 


.3007 


5.24 


124.0 


16.08 




2005, 


.3007 


5.24 


191.8 


13.38 




2005, 


.3007 


5.23 


348.9 


15.59 




2005, 


.3007 


6.07 


72.2 


14.49 




2005 


.3007 


6.61 


159.4 


12.06 




2005 


.3007 


6.72 


84.6 


14.28 




2005 


.3007 


6.89 


352.8 


13.50 




2005 


.3007 


7.11 


91.2 


8.65 




2005 


.3007 


7.38 


176.6 


16.29 




2006 


.7069 


7.42 


18.0 


15.95^ 




2005 


.3007 


7.33 


78.7 


15.73 
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Star Epoch Separation" RA* Am'' 

(arcsec) (deg) mag 

2005.3007 7.47 75.0 12.53 

2005.3007 7.71 134.8 12.62 

2005.3007 7.94 110.5 15.55 

2005.3007 8.19 198.2 15.44 

2005.3007 8.37 173.6 15.88 

2005.3007 8.60 243.7 14.96 

2005.3007 8.69 132.5 15.01 

2006.7069 9.01 318.9 15.97f 

2005.3007 9.08 141.2 12.99 

2005.3007 9.23 152.7 15.98 

2005.3007 9.32 142.4 14.95 

2005.3007 9.37 205.7 13.88 

2005.3007 9.28 337.0 12.77 

2005.3007 9.63 235.2 15.38 

2005.3007 9.56 78.8 13.78 

2006.7069 9.86 286.6 15.74^ 

2006.7069 10.03 289.1 15.75^ 

2005.3007 10.20 97.7 13.31 

2005.3007 10.18 37.3 8.81 

2005.3007 10.31 319.4 14.47 

2005.3007 10.48 308.7 15.38 

2005.3007 10.71 337.2 14.52 

2005.3007 10.92 36.3 11.98' 

2005.3007 11.36 249.0 13.24 

2005.3007 11.38 112.7 9.04 

2005.3007 11.80 40.1 14.79"= 

2006.7069 12.00 261.7 15.23f 

2005.3007 12.01 329.6 9.78 

2005.3007 12.28 153.9 14.26 

2006.7069 12.28 265.3 13.99'' 

2005.3007 12.43 114.8 9.89= 

2006.7069 12.47 71.5 14.76' 

2005.3007 12.40 100.7 13.16 

2005.3007 12.61 111.8 13.37"= 

2005.3007 13.04 16.9 14.10"= 

2006.7069 13.21 161.5 13.02' 

2006.7069 13.46 65.5 14.51' 

2005.3007 13.53 303.6 12.68"= 

2006.7069 13.59 257.1 11.76' 

2005.3007 14.12 314.4 12.55"= 

2005.3007 14.59 112.9 10.20"= 

2005.3007 14.66 17.4 12.43"= 

2005.3007 14.74 34.6 11.8P 

HD 166181 2005.2925 10.38 53.4 14.40 

2005.2925 11.21 195.8 15.04 

2005.2925 13.40 167.6 14.19^ 

2005.2925 14.46 262.8 11.93= 

HD 187748 2005.3965 5.51 325.9 15.76 

2005.3965 7.93 277.1 13.01 

2005.3965 8.02 276.7 12.18 

2005.3965 12.81 114.3 9.74 

2005.3965 13.15 321.5 12.52 

2006.7043 15.02 311.3 14.90' 

GJ 791.3 2005.3992 1.98 341.2 12.03'' 

2005.3992 2.39 51.3 13.51'' 

2005.3992 3.77 289.0 10.92'' 

2005.3992 3.80 137.6 13.99'' 

2005.3992 3.87 19.4 10.01'' 

2005.3992 4.38 201.6 14.31'' 

2005.3992 5.55 300.0 12.72'' 

2005.3992 5.96 49.6 10.28'' 

2005.3992 6.56 232.8 12.91'' 

2005.3992 6.66 155.9 13.37"' 

2005.3992 6.73 254.7 12.93"' 

2005.3992 8.01 10.3 13.03"' 

2005.3992 8.25 143.8 15.03"' 

2005.3992 8.31 71.5 15.02"' 

2005.3992 8.36 155.7 13.62"' 

2005.3992 8.89 177.3 12.86"' 

2005.3992 9.33 10.5 14.93"' 

2005.3992 9.55 276.2 15.45"' 

2005.3992 9.64 195.5 14.14"' 
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Star 


Epoch 


Separation" 


P.A.° 


Anf 






(arcsec) 


(oeg) 


mag 




2005.3992 


9.89 


347.0 


15.16'' 




2005.3992 


10.05 


255.1 


14.63'' 




2005.3992 


10.12 


201.2 


13.28'' 




2005.3992 


10.21 


310.6 


14.36'' 




2005.3992 


10.22 


328.8 


10.56'' 




2005.3992 


10.55 


166.9 


15.14'' 




2005.3992 


10.63 


80.3 


13.51'' 




2005.3992 


10.75 


326.6 


11.48'' 




2005.3992 


10.80 


57.8 


10.86'' 




2005.3992 


10.84 


51.7 


14.67'' 




2005.3992 


11.26 


243.5 


9.22'' 




2005.3992 


11.58 


315.9 


9.96'' 




2005.3992 


11.71 


14.8 


13.26'' 




2005.3992 


12.17 


18.6 


11.27'' 




2005.3992 


12.45 


46.6 


12.52'' 




2005.3992 


12.80 


274.4 


11.86'' 




2(105 3992 


13.07 


127.8 


1 1 .44'' 






13 08 


75 4 


13 58'' 




2005.3992 


13.23 


81.0 


11.26'' 




2005.3992 


13.73 


70.5 


8.91'' 




2005.3992 


14.41 


65.8 


11.04'' 




2005.3992 


15.14 


356.5 


11.98'' 




2005.3992 


15.23 


179.8 


12.02'' 


HD 201651 


2005.4867 


3.67 


201.4 


12.48 




2005.4867 


8.39 


259.1 


12.90 




2005.4867 


14.53 


331.5 


13.75= 


HD 202575 


2005.5386 


5.54 


28.5 


12.41 




2005.5386 


12.37 


168.0 


13.98 


GJ4199 


2004.6431 


9.16 


319.9 


12.24 




2004.6431 


11.76 


177.6 


10.58 


HD 206860 


2005.6069 


3.67 


60.0 


15.13 


HD 208313 


2005.4868 


2.93 


30.6 


14.78 




2005.4868 


6.24 


31.1 


9.69 




2005.4868 


9.45 


301.0 


14.41 




2005.4868 


10.45 


137.8 


16.50= 




2005.4868 


11.43 


151.6 


15.85 




2005.4868 


13.23 


121.0 


13.73 




2005.4868 


13.51 


33.2 


14.61 




2005.4868 


15.13 


63.8 


12.00 


V383 Lac 


2005.5660 


4.00 


100.0 


14.14 




2005.5660 


4.02 


79.3 


15.98 




2005.5660 


4.63 


204.6 


11.81 




2005.5660 


4.70 


207.7 


14.89 




2005.5660 


8.49 


181.6 


14.82 




2005.5660 


9.09 


110.0 


11.65 




2005.5660 


9.55 


358.6 


14.15 




2005.5660 


10.59 


93.0 


8.48 




2005.5660 


11.68 


142.3 


11.24 


HD 213845 


2005.6453 


12.85 


214.2 


14.91 


GJ 875.1 


2005.6071 


7.83 


343.9 


9.40 




2005.6071 


8.97 


260.7 


12.73 




2005.6071 


11.42 


151.3 


12.30 




2005.6071 


14.54 


15.2 


11.85= 


GJ 9809 


2006.7019 


2.10 


240.4 


15.30^ 




2005.5907 


2.18 


30.6 


12.91 




2006.7019 


3.22 


141.6 


15.55^ 




2005.5907 


3.35 


207.8 


12.31 




2005.5907 


3.55 


337.2 


11.53 




2005.5907 


3.79 


173.8 


11.92 




2006.7019 


4.51 


0.9 


16.23^ 




2005.5907 


5.37 


258.1 


14.38 




2006.7019 


5.66 


118.3 


15.35^ 




2005.5907 


6.39 


68.9 


15.47 




2005.5907 


7.01 


101.1 


9.69 




2005.5907 


7.34 


236.9 


13.49 




2005.5907 


7.69 


127.3 


11.26 




2005.5907 


7.75 


131.1 


14.21 




2006.7019 


7.92 


78.4 


15.65f 




2006.7019 


9.08 


247.3 


14.84' 




2005.5907 


9.18 


36.5 


15.03 




2005.5907 


9.23 


27.2 


14.80 




2005.5907 


9.51 


95.2 


12.35 
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Star 


Epoch 


Separation" 


RA.* 


Am'' 






(arcsec) 


(deg) 


mag 




2006.7019 


9.58 


84.1 


15.30^ 




2005.5907 


9.98 


68.5 


14.35 




2005.5907 


9.98 


121.4 


10.42 




2006.7019 


10.02 


163.4 


15.37f 




2005.5907 


10.18 


93.9 


13.69 




2006.7019 


10.73 


12.2 


15.33'' 




2005.5907 


10.64 


248.1 


13.59 




2005.5907 


10.94 


254.4 


9.90 




2005.5907 


11.23 


82.9 


13.75 




2006.7019 


11.40 


336.1 


15.18f 




2005.5907 


11.69 


155.9 


12.14 




2005.5907 


11.87 


32.2 


12.28 




2005.5907 


11.74 


291.2 


7.01 




2005.5907 


12.23 


310.2 


8.73 




2005.5907 


12.75 


66.9 


12.76 




2005.5907 


13.03 


58.0 


11.09 




2005.5907 


12.93 


332.5 


13.42 




2005.5907 


14.04 


309.7 


12.50 


HD 220140 


2005.5934 


6.14 


358.5 


15.96 




2005.5934 


15.19 


50.4 


10.62^ 


HD 221503 


2005.6646 


9.02 


234.4 


15.61'' 


GJ900 


2004.6458 


7.41 


76.0 


14.20 




2004.6458 


12.15 


150.6 


12.53 




2004.6458 


12.41 


96.4 


9.36 


GJ 907.1 


2005.6837 


7.93 


296.7 


13.68 



Note. — Target stars around which no point source was detected are omitted trom this table. Unless stated otherwise, all point sources listed were confirmed to be background 
objects using data from two epochs. 

^ Uncertainty is 0.015", see text for detail. 
^ Uncertainty is 0.2°, see text for detail. 

Uncertainties are given in Table 3, see text for detail. 

No second epoch data available. 
^ Source undetected in second epoch data. 
' Source detected in second epoch data only. 

^ Previously Imown brown dwarf companion (Potter et al. 2002; Goto et al. 2002) 

4.3. Multiple systems 

As mentioned in §2, 16 of the target stars are part of mul- 
tiple systems. As an orbital solution can potentially be de- 
termined in a reasonable amount of time for close-separation 
multiple systems, the measured properties of the systems ob- 
served that have a separation below 2" are presented in Ta- 
ble 6 as a reference for future studies. Two of the close- 
separation systems observed were resolved for the first time 
by our observations (HD 14802 and HD 166181), and a third 
system (HD 213845) was found to be a relatively large sepa- 
ration (^6") binary for which we have found no prior indica- 
tion in the literature; since HD 213845 is reported to be part 
of a binary system for the first time in this paper, its properties 
are presented in Table 6 as well. The three multiple systems 
observed for the first time in this work are discussed in more 
detail below. 

HD 14802— A source 12 ±2 times fainter than HD 14802 
was detected at a separation of 0.469" ± 0.005" and P.A. 
267.1° ±0.7° (epoch 2005.6348); the large uncertainty on 
the flux ratio is due to the peak of the primary star PSF be- 
ing saturated. Common proper motion of the pair was not 
verified but the system is likely bound given the brightness 
and close separation of the companion. The Hipparcos cata- 
log (Perryman & ESA 1997) indicates that the proper motion 
of this star is accelerating and the star is Ukely part of a binary 
system; an astrometric solution for the system was obtained 
by Gontcharov et al. (2000). The estimated period and semi- 
major axis are 25 yr and 0.5", respectively, consistent with the 



projected separation we have measured. 



HD 166181 — This star has been known for a long time to be 
a spectroscopic binary with a period of only 1.8 days (Nadal 
et al. 1974). More recently, analysis of additional radial ve- 
locity data has lead Dempsey et al. (1996) to propose that the 
system is in fact triple; a proposition which was confirmed by 
Fekel et al. (2005), who found radial velocity variations as- 
cribable to a third component with an orbit of period 5.7 year 
and eccentricity 0.765. Further, by reanalyzing Hipparcos 
data in light of this new component, these authors have found 
a new astrometric solution for the system, leading to revised 
values of parallax and proper motion (see Table 1) and to a de- 
termination of the orbital inchnation of the long-period com- 
panion. Based on their complete solution, they estimate the 
semi-major axis of the outer companion at 0.077" (2.5 AU) 
and its mass at 0.79 Mq. Our observations have resolved the 
long-period companion of this triple system. In 2005.2926, 
the companion was located at a separation of 0.065" ± 0.005" 
and PA. of 16.2° ±5.0°, and in 2006.7124, it was located at a 
separation of 0.102" ±0.003" and PA. of 51.5° ±2.0°. The 
evolution of the separation and P.A. of this source between the 
two epochs is far from that expected for an unrelated back- 
ground source and is in very good agreement with the orbital 
motion expected based on the astrometric solution of Fekel 
et al. (2005) (see Fig. 7), confirming that the source observed 
is HD 16618 IB. The flux ratio of the component Aab to com- 
ponentB is -^5.5, a contrast of -^1.85 mag. 
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TABLE 6 

Properties of close separation multiple systems 



Name 


Epoch 


Separation 


P.A. 


Brightness 


spatially resolved 






(") 


(°) 


ratio^ 


observation 



HD 14802AB 
GJ 234AB 
HD 77407 AB 
HD 102392AB 

HD 129333AB 
HD 135363AB 

HD 160934AB 

HD 166181AabB 

HD 167605AB 
HD213845AB 

GJ 900AB 

GJ 900AC 

GJ907.1AB 



2005.6348 
2005.8455 
2005.3163 
2005.3083 
2006.1929 
2005.3004 
2005.2950 
2006.1277 
2005.2953 
2006.7097 
2005.2926 
2006.7124 
2005.4017 
2005.6453 
2006.5108 
2004.6459 
2005.6864 
2004.6459 
2005.6864 
2005.6837 
2006.5411 



0.469 ± 0.005 
1.392 ±0.002 
1.702 ±0.004 
1.1 34 ±0.002 
1.1 37 ±0.002 
0.766 ± 0.002 
0.302 ± 0.002 
0.316 ±0.002 
0.213 ±0.002 
0.218 ±0.002 
0.065 ±0.005 
0.102 ±0.003 
1.182 ±0.002 
6.09 ±0.03 
6.09 ±0.03 
0.611 ±0.002 
0.673 ± 0.002 
0.733 ± 0.002 
0.722 ± 0.002 
0.789 ± 0.002 
0.775 ±0.002 



267.1 ±0.7 
44.5 ±0.1 
355.6 ±0.1 
172.5 ±0.1 
171.1±0.1 

173.0 ±0.2 
129.9 ±0.5 
131.8±0.5 

268.5 ±0.7 
271.3 ±0.7 
16.2 ±5.0 
51.5 ±2.0 
46.9 ±0.1 
129.8 ±0.4 

129.6 ±0.4 
334.3 ±0.2 
338.8 ±0.2 

344.7 ±0.2 
345.6 ±0.2 
213.5±0.2 

212.1 ±0.2 



12±2 

4.7 ±0.1 

7.4 ±0.3 

7.8 ±0.3 

14.0±0.5 
4.0 ±0.1 

2.2 ±0.1 

5.5 ±0.4 

8.2 ±0.2 
^125" 

5.7 ±0.2 

12.0 ±0.2 

1.62 ±0.05 



This work 

W. Baade (1955), c.f. Gatewood et al. (2003) 
Mugrauer et al. (2004) 
Rossiter (1955) 

Metchev & Hillenbrand (2004) 
This work'' 

Hormuth et al. (2007) 

This work 

Arribas et al. (1998) 
This work 

Martin (2003) 

Martin (2003) 

Rossiter (1955) 



^ Brightness of the primary over that of the companion, in the NIRI CH4-short filter. 

^ Independently resolved by Biller et al. (2007), evidence for co-motion is reported in this work for the first time. 
The peak of both the primary and companion is saturated in the data; the ratio quoted is an estimate based on the comparison of radial profiles. 




Fig. 7. — Verification of the physical association of the point source de- 
tected around HD 166181. Open diamonds mark the observed separation 
(top) and P.A. (bottom) of the point source at the two epochs. The solid line 
indicates the expected separation and P.A. of a distant background source as 
a function of time. The predicted separation and P.A. of HD 166181B based 
on the astrometric solution of Fekel et al. (2005) are shown as dashed lines, 
with uncertainties indicated by the shaded areas. 

HD 213845 — A bright source is visible in our data at a 
separation of 6.09" ± 0.03" and P.A. of 129.8° ± 0.4° from 
HD 213845 (epoch 2005.6453). This source did not change 
separation nor P.A. between our 2005 and 2006 observations 
(see Figure 8), indicating that it is bound to HD 213845. The 
companion is only visible in our saturated data as its separa- 
tion exceeds the field of view of the sub-array used for the 
unsaturated observations. Further, being relatively bright, the 
peak of the companion's PSF is saturated in all our data, mak- 
ing it very difficult to estimate its flux ratio to the primary 
and explaining the larger uncertainty on the separation and 
PA. quoted above. We have nevertheless estimated that the 
companion is 125 times fainter than its primary based on a 
comparison of their radial intensity profiles at radii where the 
data are in the linear regime of the detector. The companion 
was possibly detected by 2MASS, but its measured position 



Fig. 8.— Same as Figure 7 for HD 213845 



and photometry in the 2MASS point source catalog (PSC) are 
affected by confusion due to the nearby primary. Neverthe- 
less, the relative position of this source in the 2MASS PSC, 
separation of 5.55" and P.A. of 128°, is consistent with the 
star being gravitationally bound to HD 213845 as, were it 
not a bound companion, its separation should have changed 
by ^2" between the 2MASS observations and our first epoch 
observations. Although the separation of this binary system is 
well above the resolution limit of seeing-limited observations, 
we have found no prior indication of binarity in the literature. 

5. ANALYSIS AND DISCUSSION 

The detection limits determined in §4. 1 can be used to cal- 
culate an upper limit to the fraction of stars that have compan- 
ions of mass and semi-major axis inside some given intervals. 
The analysis presented in this section is largely guided by the 
work of Brandeker et al. (2006); Carson et al. (2006); Allen 
et al. (2005); and Sivia (1996). The statistical formalism for 
the analysis is presented first and various applications to our 
data are presented afterward. 
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5.1. Statistical formalism 

Consider the observation of A'^ stars enumerated by j = 
1 ...N. Let / be the fraction of stars that have at least 
one companion of mass and semi-major axis in the intervals 
['Wmin,«w] and [flmin,flmax], respectively, and pj the proba- 
bihty that such a companion around star j, if indeed it was 
there, would be detected given the detection limits of the ob- 
servations. The probabihty of detecting such a companion 
around star j is fpj, and the probability of not detecting a 
companion around this star is simply 1 -fpj- If the set {dj} 
denotes the detections made by the observations, such that dj 
equals 1 if a companion is detected around star j or else equals 
0, then the probability that the observed outcome would oc- 
cur, also called the likelihood of the data given /, is given 
by 

7=1 

According to Bayes' theorem, from the a priori probabihty 
density p{ f ), or prior distribution, and the likelihood function 
C, one may calculate p(f\{dj}), the probability density up- 
dated in light of the data, or posterior distribution: 



Pif\{dj})- 



a{dj}\f)p{f) 
j^a{dj}\f)p{f)df 



(4) 



In this study, since we have no prior knowledge about /, we 
use the most ignorant prior distribution p{f) = 1 . 

The posterior distribution p{f\{dj}) can be used to deter- 
mine a credible interval (CI) for /, bounded by /min and /max, 
for a given level of credibility a. For a case where there is no 
detection, as is the case with our survey, then clearly f^^^ = 0, 
and the upper bound of the CI is found by solving 



a ■■ 



rfm 
Jo 



p(f\{dj})df. 



(5) 



For a case where there are some detections, an equal-tail CI is 
found by solving 



1-a 



/m,„ 



p(f\{dj})df and 



1-a 



p(f\{dj})df. 

(6) 

In this work a value of a = 0.95 was chosen. 

The determination of the p/s is a critical step of this anal- 
ysis; their value depends on the detection limits of the obser- 
vations, on the ages and distances of the systems, and on the 
mass, semi-major axis, and orbital eccentricity distributions 
of the companions. In calculating the p/s it is also impor- 
tant to account properly for orbital inclination and phase as 
these affect significantly the distribution of projected separa- 
tions for an orbit of given semi-major axis. In this work, the 
p/s were calculated using a Monte Carlo approach. The mass 
and semi-major axis intervals, [mmin,?Mmax] and [amin,amax], 
were first selected. Then for each target star, 10000 planets 
were generated by sampling randomly, for each planet, the 
mass, semi-major axis, orbital eccentricity, orbital separation 
projection factor, age of the system, and underlying residual 
noise in the image. The mass and semi-major axis distribu- 
tions are left arbitrary for the moment; different possibilities 
will be explored later. For all of our calculations, the orbital 
eccentricity distribution was assumed to be that of the radial 
velocity exoplanets sample, which was approximated by a 
Gaussian function of mean 0.25, standard deviation 0.19, and 



with < e < 0.8 (Marcy et al. 2005). The orbital separation 
projection factor was sampled using the method described in 
Appendix A of Brandeker et al. (2006); this method properly 
takes into consideration orbital eccentricities, phases, and ori- 
entations. The age was sampled linearly within the range in- 
dicated in Table 1 . Given the age assigned to each planet, the 
procedure described in 4. 1 was used to convert its mass into 
a magnitude difference in the NlRl CH4-short filter. The pro- 
jected physical separation of each planet was converted into 
an angular separation based on the distance of its primary star. 
The random noise added to each planet was drawn from a 
Gaussian distribution of standard deviation equal to the mea- 
sured residual noise at the separation of the planet. The effect 
of adding this noise is to either increase or decrease the sig- 
nal that a planet would have in the residual image. The signal 
of some "a priori detectable" planets near the detection limit 
will thus decrease below the 5 a detection limit, while the sig- 
nal of some "a priori undetectable" planets will be boosted 
above the detection limit such that the appropriate detection 
completeness will result for planets of various true intensi- 
ties (see §4.1). Finally, given the sample of planets assigned 
to target j, the probability pj was calculated as the fraction 
of planets lying above the corresponding 5ct detection limits 
(c.f. Table 4). 

The above determination of the p/s yields a CI for / that 
is a function of the assumptions made on the mass and semi- 
major axis distributions. For a case where there is no detec- 
tion, it is also possible to obtain a more conservative estimate 
of /max that is vaUd for any distributions of mass and semi- 
major axis. The procedure used to do this is identical to that 
described above except for the following. Rather than popu- 
lating the whole intervals of mass and semi-major axis con- 
sidered, all planets are assigned a mass and semi-major axis 
precisely equal to mmin and flmin. respectively. Because more 
massive or more distant planets are easier to detect, the val- 
ues of p/s calculated in this manner constitute lower hmits 
to the values that would be obtained by populating the whole 
intervals assuming any specific distributions; accordingly, the 
resulting value of /max constitutes an upper limit. This ap- 
proach is perfectly legitimate as long as amax is chosen such 
that the values of p/s for any a in [amin,<3!max] are at least as 
large as those for Omin- 

It is possible to derive a simple analytic expression for /max 
for a case where there is no detection; this expression may be 
useful to estimate what the results of an ongoing survey will 
be or scale actual results for different values of or detec- 
tion probabilities. This expression may be obtained by first 
replacing each pj by the average detection probability (pj) in 
Eq. (3), and then recognizing that the likelihood function can 
be approximated by e"^/<''j) . This leads to 



/m 



-ln(l-a) 

N{Pj) 



(7) 



This approximation, valid for {p/ 3> 1, is equivalent to us- 
ing Poisson statistics rather than Binomial statistics for the 
presence of companions (c.f. Eq. 3-7 of Carson et al. 2006). 

5.2. fmaxfor arbitrary mass and semi-major axis 
distributions 

As a first analysis of the survey results, we present estimates 
of /max that are independent of the mass and semi-major axis 
distributions for Wmin =0.5, 1, 2, 3, 4, 5, 7.5, 10, and 13 Mjup, 
and for all amin between 10 and 500 AU; these estimates were 
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Fig. 9. — Mean probability of detection of a planet of given mass as a 
function of the semi-major axis of its orbit; the curves are labeled by the 
mass of the planet, in A/jup. The mean is obtained over all targets of the 
survey. 
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Fig. 10. — Upper Umits, with a credibility of 95%, on the fraction of 
stars harboring at least one companion of mass in the range 40] Mjup 
and orbit of semi-major axis in various ranges. The minimum mass, mmin, is 
indicated on each curve. For any interval, [amin , Omax] AU, of semi-major axis 
selected, the correct value of /max to read from the graph is the maximum of 
the curve within that interval. The curves shown in this graph are conservative 
upper limits that are valid for any distributions of mass and semi-major axis. 
The dotted line indicates the minimum upper limit that one could derive from 
observation of 79 stars if the probability of detection of a planet was 100% 
irrespective of its age, mass, and orbital separation. 

calculated according to the last procedure described above. 
For this analysis, and those in the next section, we have not 
considered the 6 stars with candidates for which second epoch 
observations are missing. The results obtained in this section 
are valid for any nijaax up to ~40 Mjup as no companion with 
a mass below this value was detected.^ Planet detection prob- 
abilities for each star are indicated in Table 7 for a small se- 
lection of masses and semi-major axes, while the mean planet 
detection probabilities, i.e. the average of the p/s over all 
/s, are shown in Figure 9 as a continuous function of semi- 
major axis and for a larger selection of masses. The peak 
sensitivity of the survey occurs for semi-major axes between 
50 and 200 AU; the peak values are -45% and -68% for 2 
and 5 Mjup, respectively. The survey is particularly sensitive 
to brown dwarfs (m > 13 Mjup), with a detection probabihty 
above 75% between 35 and 200 AU. A decUne in sensitiv- 
ity occurs at a separation of -^200 AU; this is consistent with 
the field of view of the observations radius) and mean 

* The previously known 40-65 Mjup binary brown dwarf companion lo- 
cated 2.6" from HD 130948 (Potter et al. 2002; Goto et al. 2002) is detected 
in our data. 
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Fig. 11. — Upper limits, with a credibiUty of 95%, on the fraction of stars 
harboring at least one planet of mass in the range [0.5, 13] Mjup, assuming 
dn/dm oc m^, and semi-major axis in various ranges. The values of (3 are 
-2 (dot-dashed line), -1.2 (solid line), and (dashed line). For any interval, 
[aminjOmax] AU, of semi-major axis selected, the correct value of /max to read 
from the graph is the maximum of the curve within that interval. The 67% 
credibility curve for /3 = -1.2 is also shown (dotted line). 

distance of the targets (22 pc). 

The results for /max are shown in Figure 10. For a semi- 
major axis interval lower bound of 50 AU, the 95% credi- 
ble planet frequency upper limits are 0.28 for 1-13 Mjup and 
50-225 AU, 0.12 for 2-13 Mj^p and 50-295 AU, and 0.057 
for 5-13 Mjup and 50-185 AU. For a semi-major axis lower 
bound of 25 AU, the upper hmits are 0.23 for 2-13 Mjup and 

25- 420 AU and 0.09 for 5-13 Mjup and 25-305 AU. For com- 
pleteness, the exercise was repeated for circular orbits and for 
a uniform distribution of eccentricity (between and 1), and 
the results obtained were very similar to those shown in Fig- 
ure 10. 

The results also indicate that no more than 0.056 of stars 
have low-mass brown dwarf companions (13 < w/Mjup < 40) 
between 25 and 250 AU. For determining the frequency of 
stars with at least one companion in the whole brown dwarf 
mass range (13 < m/Mjup < 75) over the same range of semi- 
major axis, the brown dwarf companion to HD 130948 (Pot- 
ter et al. 2002; Goto et al. 2002) must be taken into account 
explicitly. This analysis must be carried out with care as 
the semi-major axis of this companion could be significantly 
different from its measured projected physical separation of 
47 AU. It is possible to account for this uncertainty by calcu- 
lating the probability distribution of the real semi-major axis 
of the brown dwarf companion using a Monte Carlo approach 
similar to the one presented above for the calculation of the 
Pj's. Basically, the projected separation of the companion is 
fixed at s = 47 AU and its orbital eccentricity and orbital pro- 
jection factor are sampled randomly 10^ times, as described 
above. The de-projected semi-major axis is then calculated 
for each random trial and its normalized distribution over all 
trials is obtained. As the projection factor can never be larger 
than (l+Smax), where Cmax is the maximum eccentricity al- 
lowed, the semi-major axis probability distribution is equal 
to zero below ^/(l +emax); the distribution extends to infinity 
for higher values. Applied to the current case, this calculation 
indicates that at a 95% credible interval for the semi-major 
axis of the binary brown dwarf companion to HD 130948 is 

26- 157 AU. We thus posit that our observations have resulted 
in one detection in the semi-major axis interval 25-200 AU 
and mass interval 13-75 Mjup; then using the procedure de- 
scribed in the previous section and Eq. (6), the 95% credi- 
ble interval for the frequency of stars with at least one brown 
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Fig. 12. — Upper limits, with a credibility of 95% (top panel) or 67% (bottom panel), on the fraction of stars harboring at least one giant planet of mass in 
the range [0.5, 13] Mjup, assuming dn/dm oc m~'-^, and orbit of semi-major axis in the range [amm,ama}!] AU, assuming dn/da oc a''. The abscissa indicates the 
upper bounds (omax) of the semi-major axis intervals, while the lower bounds (Omin) are 10 AU (solid lines), 25 AU (dotted lines), and 50 AU (dashed lines). The 
top, middle, and bottom curves in each set of three curves are for 7 = -1, 0, and 1, respectively. 



dwarf companion in the semi-major axis interval 25-250 AU 
is 0.019^g Qj5. This result is consistent with the upper limit of 
0. 12 (95% credibility) reported by Carson et al. (2006) for the 
25-100 AU semi-major axis interval and also with the frac- 
tion of 0.068!° (95% credibihty) reported by Metchev & 
Hillenbrand (2004) for the range 30-1600 AU. For smaller 
semi-major axes, our results indicate that, with a credibility 
of 95%, the fraction of stars with at least one brown dwarf 
companion in the range 10-25 AU is less than 0.20, and less 
than 0.10 for the range 15-25 AU. 

5.3. /max for specific mass and semi-major axis distributions 

In this section we derive first an upper limit to the fraction 
of stars harboring at least one planet in the single mass inter- 
val [0.5, 13] Mjup, assuming that the mass distribution follows 
dn/dm oc m"'^. The mass distribution adopted is based on a 
statistical analysis of the RV results that properly accounts for 
the detection sensitivity reached for each star (A. Gumming et 
al. 2007, in preparation) and is formally valid only for planets 
with semi-major axis below ^3 AU; here it is blindly extrap- 
olated to larger semi-major axes. For comparison, a simple fit 
of the mass distribution of the RV exoplanets sample yields 
dn/dm oc m~'-^ (Butler et al. 2006). For this calculation the 
whole mass interval is populated according to the distribution 
stated, but all planets are assigned a value Omin for the semi- 
major axis, so as to make the results independent of its dis- 
tribution. The calculation was made for all amin between 10 
and 500 AU. The results are shown in Figure 1 1 . With a cred- 
ibility of 95%, the fraction of stars having at least one planet 
of mass in the range [0.5,13] Mjup and semi-major axis in 



[10,500], [25,340], and [50,230] AU is less than 0.57, 0.17, 
and 0.10, respectively. For reference, results of the same anal- 
ysis assuming dn/dm oc m'^, with (3 = and -2, are presented 
also in Figure 1 1 . As expected, a smaller /3 leads to larger val- 
ues of /max because a larger fraction of planets have a smaller 
mass, while a larger value of /3 has the opposite effect. 

Next we calculate upper limits for the same mass interval 
by assuming further that the distribution of semi-major axes 
follows 6n/6a oc a'^, for 7 = -1, 0, and 1. This range of 
power-law index includes the value of 7 = -0.67 found by A. 
Gumming et al. (2007, in preparation) for the RV exoplanets 
sample within the range 0.03-3 AU. We have done the cal- 
culations for amin=10, 25, and 50 AU, and for all Omax in the 
range [amin+5,500] AU; the results are presented in Figure 12 
and the corresponding planet detection probabiUties for each 
star are shown in Table 7 for a selection of semi-major axis 
intervals. Figure 13 illustrates how a synthetic population of 
planets based on the above distributions compares with our 
detection limits. For 7 = -1, the 95% credible upper limits 
to the fraction of stars with at least one planet of mass in the 
range [0.5, 13] Mjup are 0.28 for the semi-major axis range 
10-25 AU, 0.18 for 10-50 AU, 0.13 for 25-50 AU, 0.11 for 
25-100 AU, and 0.093 for 50-250 AU. Slightly smaller val- 
ues of /max are found for larger values of 7, as such indices 
would place more planets at larger separations where they 
would have been easier to detect with our observations. For 
the larger values of flmin> the value of 7 has very little effect 
on the upper limit found as, irrespective of the value of 7, 
the majority of planets are located at separations for which 
the sensitivity of the observations is high. Overall, the weak 
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Fig. 13. — Detection limits (5u, solid line) and synthetic population of 
planets (dots) for the star HD 166181. A planet mass distribution following 
dn/dm oc m~'^ inside 0.5-13 Mj^p and a semi-major axis distribution follow- 
ing dn/dacx a"' inside 10-300 AU were used. For this particular example, 
the planet detection probability pj is 30%. 

dependence of /max on 7 implies that the semi-major axis dis- 
tribution (i.e. 7) cannot be constrained from our results. 



As one may worry that the population of planets around M 
dwarfs is different from that around earlier-type stars, because 
of smaller disk masses for example, we derive an estimate of 
/max by excluding the M dwarfs from the statistical analysis. 
This estimate is obtained using Eq. (7) and the values of the 
last three columns of Table 7; it is thus valid for /3 = -1 .2 and 
7 = -1. Excluding M dwarfs from the sample leaves 64 stars 
and results in average detection probabilities (pj) of 0.070, 
0.229, and 0.385 for 10-25 AU, 25-50 AU, and 50-250 AU, 
respectively. The corresponding 95% credible upper limits to 
the fraction of stars with planets are then 0.67, 0.20, and 0.12. 
The effect is quite significant at the smallest orbilal separa- 
tions, where M dwarfs provide good sensitivities due to their 
smaller luminosity, smaller average distance, and younger av- 
erage age. Similarly, the population of planets in stellar mul- 
tiple systems may be different from that in single systems. 
Excluding multiples from the sample also leaves 64 stars and 
yields values of (pj) of 0.139, 0.299, and 0.389 for 10-25 AU, 
25-50 AU, and 50-250 AU. The corresponding upper limits 
to the fraction of stars with planets are 0.34, 0.16, and 0.12; 
the effect is thus rather small in this case. 



TABLE? 
Planet detection probability 



Name 






2Mj„p 










5Mj„p 








0.5-13 Mj„p, 


/3 = -1.2,7 = 


-1 


25 AU 


50 AU 


100 AU 


25 AU 


50 AU 


100 AU 


10- 


25 AU 


25- 


50 AU 


50-250 AU 


HD 166 





002 





261 





748 





693 





936 





985 





120 





339 





448 


HD691 





000 





078 





403 





113 





730 





944 





022 





191 





452 


HD 1405 





009 





602 





919 





540 





906 





979 





076 





348 





654 


HD 5996 





000 





035 





164 





134 





575 





900 





030 





176 





344 


HD 9540 





001 





041 





117 





112 





370 





641 





028 





148 





269 


HD 10008 





000 





010 





282 





161 





792 





957 





044 





212 





419 


GJ82 





921 





983 





963 





978 





993 





967 





575 





788 





648 


HD 14802 





000 





000 





000 





000 





000 





000 







n 
u 


\JV\J 





000 


HD 16765 





000 





019 





231 





149 





730 





944 





029 





192 





389 


HD 17190 





001 





016 





041 





037 





115 





199 





008 





059 





151 


HD 17382 





014 





138 





322 





296 





752 





949 





064 





236 





413 


HD 17925 





867 





972 





843 





947 





987 





848 





454 





697 





508 


HD 18803 





000 





000 





000 





000 





000 





008 





001 





019 





087 


HD 19994 





000 





000 





000 





000 





000 





000 





000 





005 





062 


HD 20367 





000 





003 





222 





012 





486 





894 





009 





123 





393 


2E759 





001 





197 





648 





396 





872 





973 





089 





282 





490 


HD 22049 





000 





000 





000 





935 





338 





068 





273 





257 





035 


HIP 17695 





772 





951 





988 





944 





987 





996 





386 





614 





633 


HD 25457 





000 





030 





635 





148 





797 





956 





022 





203 





449 


HD 283750 





369 





864 





971 





840 





966 





991 





202 





506 





656 


GJ 182 





530 





906 





979 





828 





964 





991 





222 





610 





873 


GJ281 





021 





310 





548 





789 





956 





989 





169 





375 





416 


GJ285 





982 





992 





405 





987 





993 





407 





811 





961 





425 


HD 72905 





000 





000 





052 





157 





804 





955 





046 





209 





322 


HD 75332 





000 





026 





472 





029 





624 





924 





Oil 





148 





456 


HD 77407 





007 





381 





854 





282 





828 





965 





027 





266 





649 


HD 78141 





116 





719 





942 





708 





938 





986 





121 





413 





601 


HD 82558 





373 





867 





971 





825 





962 





990 





184 





505 





661 


GJ393 





971 





992 





511 





981 





994 





513 





664 





789 





402 


HD 90905 





000 





060 





553 





073 





716 





940 





014 





173 





470 


HD 91901 





000 





000 





005 





001 





020 





058 





001 





009 





045 


HD 92945 





003 





444 





886 





517 





902 





978 





079 





322 





530 


HD 93528 





000 





020 





371 





014 





570 





915 





009 





132 





431 


GJ402 





958 





942 





271 





991 





951 





274 





581 





672 





225 


HD 96064 





064 





642 





927 





664 





929 





984 





100 





383 





598 


HD 97334 





007 





252 





610 





513 





902 





978 





063 





300 





489 


HD 102195 





000 





002 





015 





010 





057 





116 





003 





027 





097 


HD 102392 





000 





004 





018 





014 





064 





124 





005 





030 





100 


HD 105631 





000 





000 





000 





000 





000 





000 





000 





017 





107 


HD 107146 





000 





197 





798 





286 





845 





966 





035 





250 





533 


HD 108767B 





001 





122 





456 





169 





774 





953 





038 





214 





468 


HD 109085 





000 





000 





000 





000 





000 





093 





001 





022 





128 


BD+60 1417 





749 





946 





987 





938 





985 





996 





384 





612 





645 


HD 111395 





000 





000 





000 





002 





175 





532 





018 





118 





207 


HD 113449 





005 





488 





894 





604 





920 





982 





096 





340 





548 


GJ 507.1 





000 





010 





023 





033 





093 





141 





Oil 





053 





100 


HD 116956 





000 





051 





237 





174 





715 





942 





043 





200 





375 
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Name 




2Mjup 






5Mjup 




0.5-13 Mj„p, /3 = - 


-1.2, 7 = -1 


— — 7m — 

25 AU 


— CM Air — 

50 AU 


— 1 JUk All — 

100 AU 


— — ATT — 

25 AU 


— r/v A 1 f — 

50 AU 


— M\l\ All — 

100 AU 


10-25 AU 


25-50 AU 50-250 AU 


HD 118100 


0.155 


0.802 


0.956 


0.810 


0.959 


0.989 


0.191 


0.452 


0.584 


GJ 524.1 


0.001 


0.012 


0.025 


0.034 


0.100 


0.147 


0.01 1 


0.058 


0.116 


HD 124106 


0.000 


0.000 


0.000 


0.000 


0.002 


0.105 


0.003 


0.056 


0.163 


HD 125161B 


0.000 


0.008 


0.022 


0.016 


0.064 


0.116 


0.004 


0.031 


0.092 


HD 129333 


0.014 


0.296 


0.799 


0.246 


0.817 


0.962 


0.028 


0.254 


0.582 


HD 130004 


0.001 


0.011 


0.026 


0.037 


0.099 


0.154 


0.012 


0.056 


0.114 


HD 130322 


0.000 


0.000 


0.000 


0.000 


0.000 


0.005 


0.001 


0.016 


0.104 


HD 130948 


0.023 


0.451 


0.876 


0.584 


0.916 


0.981 


0.071 


0.336 


0.567 


HD 135363 


0.005 


0.370 


0.858 


0.333 


0.856 


0.969 


0.045 


0.284 


0.568 


HD 141272 


0.001 


0.167 


0.512 


0.573 


0.915 


0.980 


0.086 


0.304 


0.462 


HD 147379B 


0.168 


0.264 


0.252 


0.815 


0.961 


0.872 


0.241 


0.372 


0.279 


GJ628 


0.077 


0.056 


0.013 


0.468 


0.333 


0.067 


0.209 


0.240 


0.048 


HIP 81084 


0.002 


0.364 


0.863 


0.494 


0.898 


0.977 


0.085 


0.315 


0.528 


HD 160934 


0.500 


0.900 


0.978 


0.855 


0.968 


0.992 


0.231 


0.574 


0.780 


HD 162283 


0.002 


0.019 


0.035 


0.049 


0.125 


0.200 


0.016 


0.072 


0.140 


HD 166181 


0.000 


0.058 


0.508 


0.119 


0.756 


0.949 


0.021 


0.193 


0.468 


HD 167605 


0.000 


0.003 


0.015 


0.004 


0.043 


0.093 


0.002 


0.018 


0.078 


HD 187748 


0.000 


0.126 


0.689 


0.181 


0.799 


0.958 


0.036 


0.223 


0.516 


HD 197481 


0.910 


0.980 


0.808 


0.913 


0.980 


0.810 


0.538 


0.897 


0.622 


HD 201651 


0.000 


0.004 


0.015 


0.008 


0.043 


0.089 


0.002 


0.019 


0.073 


HD 202575 


0.009 


0.102 


0.204 


0.262 


0.651 


0.921 


0.070 


0.218 


0.333 


GJ 4199 


0.518 


0.904 


0.978 


0.841 


0.965 


0.991 


U.zi J 


0.567 


0.799 


HD 206860 


0.000 


0.018 


0.255 


0.427 


0.882 


0.974 


089 


0.262 


0.413 


HD 208313 


0.012 


0.106 


0.200 


0.296 


0.676 


0.923 


0.082 


0.228 


0.351 


V383 LAC 


0.024 


0.443 


0.869 


0.584 


0.915 


0.980 


0.105 


0.347 


0.567 


HD 213845 


0.000 


0.001 


0.172 


0.056 


0.710 


0.938 


0.016 


0.168 


0.394 


GJ 875.1 


0.214 


0.793 


0.954 


0.917 


0.982 


0.994 


0.284 


0.463 


0.479 


GJ876 


0.059 


0.054 


0.013 


0.372 


0.333 


0.068 


0.160 


0.219 


0.053 


GJ 9809 


0.811 


0.960 


0.990 


0.920 


0.982 


0.994 


0.396 


0.717 


0.893 


HD 220140 


0.150 


0.736 


0.945 


0.771 


0.952 


0.988 


0.136 


0.426 


0.573 


GJ900 


0.008 


0.443 


0.874 


0.644 


0.928 


0.984 


0.116 


0.353 


0.491 


GJ 907.1 


0.000 


0.000 


0.007 


0.009 


0.528 


0.906 


0.015 


0.134 


0.344 


Mean 


0.144 


0.304 


0.456 


0.388 


0.635 


0.696 


0.122 


0.284 


0.392 



6. SUMMARY AND CONCLUSION 

In this paper, we have presented the results of the Gemini 
Deep Planet Survey, a near-infrared adaptive optics search for 
giant planets on orbits of 10-300 AU around nearby young 
stars. The use of angular differential imaging at the Gemini 
North telescope has enabled us to reach the best sensitivities 
to date for detecting giant exoplanets with projected separa- 
tions above ~0.7". The typical detection limits (5c7) reached 
by the survey, in magnitude difference between an off-axis 
point source and the central star, are 9.5 at 0.5", 12.9 at 1", 
15 at 2", and 16.5 at 5", sufficient to detect planets more mas- 
sive than 2 Mjup with a projected separation of 40-200 AU 
around a typical target star. More than 300 faint point sources 
have been detected around 54 of the 85 stars observed, but 
observations at a second epoch have revealed changes in sep- 
aration and P.A. of these point sources relative to the target 
stars that are all consistent with those expected from unre- 
lated background objects. The observations made as part of 
this survey have resolved the stars HD 14802, HD 166181, 
and HD 213845 into binaries for the first time. 

We have presented a statistical analysis of the survey results 
to derive upper limits to the fraction of stars having planetary 
companions. This analysis indicates that the 95% credible up- 
per limit to the fraction of stars harboring at least one planet 
more massive than 2 Mjup with an orbit of semi-major axis in 
the range 25-420 AU or 50-295 AU is 0.23 or 0.12, respec- 
tively, independently of the mass and semi-major axis distri- 
butions of the planets; for planets more massive than 5 Mjup, 
the upper Umits are 0.09 for 25-305 AU and 0.057 for 50- 
185 AU. It was also found that less than 0.056 of stars have 
low-mass brown dwarf companions (13 < m/Mjup < 40) be- 
tween 25 and 250 AU (see Figure 10); this upper limit is set 



by the sample size only as the sensitivity of the observations 
to brown dwarfs is very good. Considering the whole brown 
dwarf mass range, the 95% credible interval for the frequency 
of stars with at least one brown dwarf companion in the semi- 
major axis interval 25-250 AU is 0.019![5;[|^^. Assuming a 
mass distribution following dn/dm oc m"'-^, the results indi- 
cate that with a credibility of 95% the fraction of stars having 
at least one planet of mass in the range 0.5-13 Mjup and semi- 
major axis in the range 25-325 AU is less than 0.17, and less 
than 0.10 for the range 50-220 AU. Assuming further a semi- 
major axis distribution following dn /da oc a"' , the upper Um- 
its to the fraction of stars with planets are 0.28 for the range 
10-25 AU, 0.13 for 25-50 AU, and 0.093 for 50-250 AU. 

The work presented in this paper constitutes a first step to- 
ward the detection of the population of "outer" giant planets 
around other stars. Such a study, which is complementary to 
RV searches in terms of orbital separation, is necessary to im- 
prove our understanding of the various mechanisms that could 
generate planets on orbits of tens to hundreds of AU, such as 
in situ formation triggered by collisions of stars with proto- 
planetary disks or orbital migration induced by gravitational 
scattering in multiple planet systems. While the upper limits 
we have found rule out an important increase in the population 
of planets at large separations compared to the known popu- 
lation of planets below 3 AU, our sample size and the sensi- 
tivities we haved reached are insufficient to tell if the above 
mechanisms operate at all, and a fortiori which one is domi- 
nant. Future searches reaching better sensitivities and targeted 
at a larger sample of stars will be necessary to answer these 
questions. 

Considerable efforts are currently deployed by major ob- 
servatories to develop instruments dedicated to the search of 
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giant exoplanets around nearby stars. The Gemini Planet Im- 
ager (GPI, Gemini Telescope, Macintosh et al. 2006) and 
the Spectro-Polarimetric High-contrast Exoplanet Research 
instrument (SPHERE, Very Large Telescope, Dohlen et al. 
2006) are good examples; they should see their first light in 
around 2010. These complex instruments will ally an extreme 
AO system to correct atmospheric wavefront errors to un- 
precedented levels of accuracy, a calibration system to correct 
instrumental quasi-static aberrations, a coronagraph to sup- 
press the coherent on-axis stellar light, and differential imag- 
ing capabilities enabled by either multi-channel cameras or in- 
tegral field spectrographs. The expected performance of these 
instruments, e.g. a contrast better than 17.5 mag at a sepa- 
ration of 0.5" for GPI (Macintosh et al. 2006), should allow 
detection of planets of 1 Mjup aged less than 100-200 Myr at 
separations of 5-50 AU, significantly improving on the work 
presented here. These efforts should uncover the population 
of outer giant planets, if they exist, or place sufficient con- 
striiints on their existence to rule out the mechanisms that 
could generate them. In less than a decade the James Webb 
Space Telescope will allow similar studies to be done for rela- 
tively nearby M-type primaries, which are too faint for operat- 
ing the wavefront sensor of extreme adaptive optics systems. 
Given all of the projects that should unfold in the next few 
years, the coming decade promises to be extremely exciting 
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